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Abstract 

This review is devoted to generalization of dynamical mean-field theory (DMFT) 
for strongly correlated electronic systems towards the account of different types of 
additional interactions, necessary for correct physical description of many experi- 
mentally observed phenomena in such systems. As additional interactions we con- 
sider: (1) interaction of electrons with antiferromagnetic (or charge) fluctuations of 
order parameter in high-Tc superconductors leading to the formation of pseudogap 
state, (2) scattering of electrons on static disorder and its role in general picture of 
Anderson-Hubbard metal-insulator transition, (3) electron-phonon interaction and 
corresponding anomalies of electronic spectra in strongly correlated systems. Pro- 
posed DMFT+S approach is based on taking into account above mentioned inter- 
actions by introducing additional self-energy E (in general momentum dependent) 
into conventional DMFT scheme and calculated in a self-consistent way within the 
standard set of DMFT equations Here we formulate general scheme of calculation of 
both one-particle (spectral functions and densities of states) and two-particle (op- 
tical conductivity) properties. We examine the problem of pseudogap formation, 
including the Fermi arc formation and partial destruction of the Fermi surface, 
metal-insulator transition in disordered Anderson-Hubbard model, and general pic- 
ture of kink formation within electronic spectra in strongly correlated systems. 

DMFT+S approach is generalized to describe realistic materials with strong 
electron-electron correlations based on LDA+DMFT method. General scheme of 
LDA+DMFT method is presented together with some of its applications to real sys- 
tems. The LDA+DMFT+S approach is employed to modelling of pseudogap state 
of electron and hole doped high-T c cuprates. Comparison with variety of ARPES 
experiments is given. 

PACS: 71.10.Fd, 71.10.Hf, 71.20.-b, 71.27.+a, 71.30.+h, 72.15.Rn, 74.72.-h 
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1 Introduction 



Strongly correlated electronic systems (SCS), which are mainly realized in a number com- 
pounds of transition or rare-earth elements with partially filled 3d, 4/ and 5/ shells, for 
more than half a century attract a lot of interest of scientists because of their unusual phys- 
ical properties and difficulties in their theoretical description. Problem of metal-insulator 
phase transition, observed in many transition metal oxides, heavy fermions systems, with 
great variety of different phase transitions and related phenomena, manganites with giant 
magnetoresistance - all these systems become a subject of great attention for both exper- 
imentalists and theorists. Perhaps the most significant development in this area was the 
discovery of high temperature superconductivity in copper oxides, which provoked a new 
wave of interest in the synthesis and description of such systems. 

As already stressed above, the diversity of physical phenomena in all these compounds 
is due to partially filled 3d, 4/ and 5/ shells. Strong interaction of electrons within narrow 
bands belonging to these orbitals shells with each other or with itinerant electrons of outer 
shells is basically responsible for unique properties of these systems. Early qualitative 
ideas formulated by Mott [1] were further developed in many theoretical works to follow, 
forming a new area of modern theory of condensed matter. There are now thousands of 
such papers and many new reviews and books are regularly published [2]. 

Of course, nowadays, a single review can not cover all this area of research and the aim 
of authors is rather modest. The object of this work is description of a number of theoret- 
ical approaches formulated during recent years to account for some additional interactions 
which are important for the coherent description of strongly correlated systems. Here we 
speak not only about "external" perturbations due to interaction of correlated electrons 
with bosonic excitations such as phonons, spin fluctuations or scattering of electrons by 
disorder, but also about attempts to improve most developed and widely used theoretical 
approaches such as dynamical mean- field theory (DM FT) [3j HJ G3 E] . 

These tasks are closely related with variety of topical problems under active study at 
present, such as the nature of the pseudogap state of high-T c cuprates, the evolution of 
their Fermi surfaces upon doping with formation of so called "Fermi arcs", observed in 
ARPES experiments, the problem of formation of kinks in electronic spectrum, the general 
problem of metal-insulator transition in disordered systems, due to mutual interference 
of strong correlations and Anderson localization. In this review to some extent we deal 
with all of these problems. 

During last years the general theory of strongly correlated systems based on DMFT 
practically merged with the so called "first-principle" approaches to calculations of elec- 
tronic spectra of real solids (LDA+DMFT approach), with significant successes already 
achieved [21 El- To this end in our review we devote some attention to first attempts 
of generalizing these approaches towards the account of the above mentioned physical 
effects. 
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2 Strongly correlated systems and dynamical mean- 
field theory (DMFT). 

2.1 Hubbard model and basics of DMFT. 

Starting with pioneering works of Hubbard [8] in the early 60th the simplest model to 
describe strongly correlated systems is the so called Hubbard model. One band Hubbard 
model Hamiltonian is: 

H = -t °L c ja + u Yl n ^ n ^ (!) 

(ij)a i 

where t > - nearest neighbors hopping amplitude, U - single site repulsion, n ia = c\ a c ia . - 
particle number operator on site i, Ci a (c\ a ) - annihilation (creation) electron operators for 
spin a. The model has only two competing energy parameters. Parameter t defines kinetic 
energy of electron and facilitates intersite hoppings (derealization), while parameter U 
defines potential energy and favors localization of an electron on a lattice cite. Energy 
bands formed by 3d, 4/ and 5/ orbitals are rather narrow, thus quite often kinetic and 
potential energy are of the same order of magnitude (t ~ U). In this case there is no small 
parameter in the model and it is impossible to build any kind of perturbation theory. This 
fact alone leads to all the difficulties in description of SCS even for such oversimplified 
model. 

Almost for 30 years there were no satisfactory approaches to analyze SCS. It seemed 
that theory of these systems will forever remain fragmentary and semiquantitative. The 
breakthrough came in 1989 in the work by Metzner and Vollhardt [9J. They suggested 
formal consideration of the system of interacting electrons in large space dimensions d — > 
oo (or in a lattice with large coordination numbers — > oo 0). 

Employing this limit it is possible to neglect spatial fluctuations in the systems while 
full local dynamics is preserved. In Ref. [9] it was shown that in the limit of infinite 
spatial dimensions (or more precise infinite coordination number) the main role is played 
only by local contributions to self-energy part of full interacting Green's function. All non 
local contributions are proportional to 1/ y/z ~ 1/ yd and can be dropped. In this limit 
electron self-energy does not depend on momentum k, and is a function of frequency only 
(real (to) or Matsubara (co n ) one) d: 

E CT (k,u;) = £ CT (u;). (2) 

This statement is the main simplification appearing in the limit of (i — )■ oo [21 El El [7] . 

In Fig. [TJ we show the "skeleton" diagrams of DMFT self-energy S. Wavy lines 
represent local (Hubbard) interaction U, full lines represent local Green functions Gu. 
Strictly speaking in the limit of d — > oo the self-energy is not only a local one, but also 
in each vertex of "skeleton" diagram, only one particular site enters, e.g. the i-th ones, 
as shown in Fig. [TJ Thus this self-energy is a functional Sjj iCr = F[Gu ja ] of interacting 
local Green's function Gu t(7 . Generally it is not sufficient to make our problem the fully 

^or hypercubic lattice z — 2d and these two limits practically coincide. However, even for three 
dimensional lattices z could be quite large, for example, in body centered lattice z = 8 and for face 
centered z = 12. To this end it is more correct to speak about limit of large z. 

2 Large coordination numbers allow one to apply this approximation rather successfully even for rather 
small d 
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Figure 1: "Skeleton diagrams of local self-energy E in DMFT. 
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Figure 2: Within the DMFT lattice Hubbard model maps to interacting electrons on a single 
site ("impurity"), surrounded by the fermionic bath defining dynamical (Weiss) mean- field G(lj). 

local one, since interacting Green function Gi^ a is still nonlocal. Then a question arises, 
whether it is possible to choose purely local non interacting (in the absence of U) problem 
with completely equivalent self-energy? Surely it can be done! Let G a (u>) to be "bare" 
dynamical Green function of such local problem without Coulomb interaction U, while 
Gdaipj) and S rf(7 (a;) are corresponding interacting Green's function and self-energy. If one 
can guarantee the equality G^oj) = Ga^(u), then corresponding self-energies are also 
equivalent, because the structure of diagrams of weak coupling U perturbation theory 
is totally preserved, which means that self-energy of the local problem is defined by the 
same functional S d(T = F[Gda]- But Q a , Gd a and S rf(7 are connected through the Dyson 
equation, which immediately gives us the "bare" dynamical Green function of the local 
problem. 

^VH^H + ^M; (3) 

Thus the lattice Hubbard model in the limit of d — > oo exactly maps onto purely 
local dynamical problem. Physically it corresponds (as shown in Fig. [2]) to the problem 
of interacting electrons on a single "Anderson impurity" in a "bath" and interaction 
with the bath is contained in dynamical mean-field Quite often, in analogy with 

molecular field theory in magnetism, this field is called "Weiss field". This explains the 
name of such an approach - dynamical mean-field theory (DMFT). 

This purely dynamic problem is still quite complicated. However, the problem is 
equivalent to a single impurity Anderson model (SIAM) [10]. This model can be studied 
in detail by a number of different methods and its physics is now well understood. For 
this model there are well developed approximate analytical methods such as iterative 
perturbation theory (IPT) [3] and non-crossing approximation (NCA) (TTJ [T2| , but most 
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Figure 3: DMFT(NRG) densities of states at half-filling obtained for semielliptic "bare" density 
of states for different values of U. 



remarkable is the possibility to solve this model by exact numerical methods like quantum 
Monte-Carlo (QMC) [13] or numerical renormalization group (NRG) [T4l [T5] . Solution of 
an effective SIAM employing any of these methods, usually called an "impurity solver", 
completes the general scheme of DMFT approach. 

Apparently, today DMFT is the most elaborate and reliable theoretical method to 
describe SCS. In its framework the so called three-peak structure of the density of states 
of SCS was obtained for the first time [5], consisting of the central (quasiparticle) peak on 
the Fermi level and two wide maxima, corresponding to upper and lower Hubbard bands. 
Also the reliable theoretical description of Mott-Hubbard metal-insulator transition was 
obtained. In Fig. [3] we show DMFT(NRG) densities of states of the half-filled Hubbard 
model with semielliptic "bare" density of states with bandwidth 2D. As correlation 
strength U grows the density of states demonstrates the formation of characteristic three- 
peak structure and further increase of U leads to a collapse of the quasiparticle peak at 
U/2D 1.5, leading to metal-insulator transition. 

It turns out that within DMFT it is also possible to investigate some two-particle 
properties. In particular, it is quite easy to obtain dynamic optical conductivity (5J H]. 
During the recent years DMFT approach was generalized to describe realistic SCS merging 
it with Ll ab initio" one-electron density functional theory in local density approximation 
(DFT/LDA), leading to the combined computational scheme of LDA+DMFT [Tol HTf [T5 1 
Q3H [2D] , which will be described later. 

Despite all the obvious advantages of DMFT this approach has a number of shortcom- 
ings. Namely, as we stressed above, all non-local correlations are completely neglected. A 
number of cluster generalizations of DMFT [211 122] were proposed recently to overcome 
this drawback. However, all these methods are quite computer time consuming and are 
rather restricted in a cluster size and with respect to their generalization to multi-orbital 
case. Also in these approaches it is quite difficult to investigate two-particle properties. 
To overcome these difficulties we proposed recently [231 [2U [25] the new generalization of 
the conventional DMFT, allowing to consider non-local correlations or additional (with 
respect to the Hubbard one) interactions (in principle of any kind), remaining within a 
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single impurity DFMT picture and preserving self-consistent set of DMFT equations. 



2.2 Generalized DMFT+E approach. 

The main idea of the new approach is to use the exact in the limit of d — > oo DMFT 
solution as a "high energy" zeroth order approximation, describing electronic spectra on 
a large energy scale of the order of bandwidth or U value, while low energy scale details 
caused by non-local effects or by effects of interaction of correlated electrons with different 
collective modes (e.g. phonons or spin fluctuations) are to be taken into account within 
some kind of perturbation theory, conserving, as far as possible, the general structure of 
DMFT equations. Actually, such a scheme can be realized in rather simple way [2 ^ 1 ^ 125]. 

To be more specific, in the following we consider the standard one band Hubbard 
model. Generalizations towards multi orbital case are also possible. Main assumption of 
our approach is to choose Matsubara lattice Green's function as: 

Gk(iw) = -— 1 , u = nT(2n + 1), (4) 

iu + (j, — £[K) — — Eic(zo;J 

where \x is the chemical potential, H(iou) - local DMFT self-energy due to Hubbard in- 
teraction and £k(iu;) - some "external" (generally non-local, momentum dependent) self- 
energy. This last contribution can arise from interaction of correlated electrons with some 
"additional" collective modes or order parameter fluctuations appearing in the Hubbard 
model itself, or from any other interactions (fluctuations) external with respect to the 
standard Hubbard model. For example these can be phonons or scattering by impurities, 
when it is actually local (momentum independent). 

One should emphasize that Ek(icu) can contain local (momentum independent) con- 
tribution even if the self-energy is considered in the framework of the Hubbard model. 
However this contribution disappears in the infinite spatial dimensions limit of d — > oo 
and is not accounted within the conventional DMFT, so that within our approach we not 
encounter any double counting problem even in this case. This question does not come 
out at all for self-energy E k (iw) caused by "external" interactions. 

More important is that our assumption of the additive form of a self-energy Yt(iu) + 
Ek(io') implicitly corresponds to the neglect of interference of the local (DMFT) and non- 
local contributions. In Fig. 0] typical "skeleton" diagrams for self-energy of DMFT+E 
approach are given. First two terms are local DFMT self-energy diagrams, two diagrams 
in the middle show contributions to non-local part of self-energy from "additional" inter- 
actions with collective modes or order parameter fluctuations, and the last diagram (b) 
is an example of diagram with interference between local and non-local parts which is 
neglected. Indeed, once we neglect such interference (i.e. diagram shown in Fig. H^b)) 
the total self-energy is defined as a simple sum of these two contributions shown in Fig. 
HJ Two last diagrams in Fig. H^a) are an example of "skeleton" diagrams for non-local 
self-energy, where full line is the Green's function (j3J) and dashed line corresponds to 
an "additional" interaction with collective modes or order parameter fluctuations. 

Finally, diagrammatic structure of the local self-energy remains identical to that of the 
standard DMFT and we obtain the following self-consistent equations of the generalized 
DMFT+E approach 



1. Start from some initial guess for the local self-energy E(iu;), e.g., E(zcu) = 0. 
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Figure 4: Typical "skeleton" self-energy diagrams of DMFT+S approach. 

2. Calculate self-energy T,^(ioj) in the framework of some (approximate) scheme, tak- 
ing into account interaction of correlated electrons with collective modes or order 
parameter fluctuations, which, in general, can depend on Y}(iu) and [i. 

3. Calculates local Green's function: 

1 x 1 

Gll{lUj) = N S icj + fj,-e(k)-Z(iu)-Z k (iuy (5) 

4. Define the "Weiss field" as: 

G^(iu) = E(iuj)+G^(iu;). (6) 

5. Using some "impurity solver" calculate the single-particle Green's function of an 
effective single impurity Anderson model, i.e. compute the following integral over 
Grassmann variables cf c and c^: 

G d (r - T ') = — I D4 a Dc ia c ia {T)4 a (j') exp(-S cS ), (7) 

^eflf J 

where an effective action for the fixed site ( "Anderson impurity" ) % is: 

S e s=- drx / dr 2 c ia (Ti)Qo 1 (T 1 - t 2 )c+ (t 2 ) + / rfrf/n it (r)n 4 (r) , (8) 
Jo Jo Jo 

with "partition function" Z cS = J Dcf a Dci a exp(—Seg), and (3 = T _1 . 

6. Define new local self-energy as: 

sh^^h-g^h. (9) 

7. Using this self-energy as an "initial" on step 1, continue the loop procedure until it 
converges to 

G u (iu>) = G d (iu). (10) 

with a given accuracy. 

At the end we obtain the final Green's function in the form of Eq. (j3J), where H(iu) 
and Sk(^w) are self-energies coming out of our iterative procedure. 

Success of such approach (as well as its main drawback) is connected with an additive 
form of total self-energy (neglect of interference between different contributions) in Eq. 
(j3J). This allows one to preserve self-consistent set of equations of the standard DMFT. 
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However there are two significant distinctions from conventional DMFT. First of all, the 
local Green's function of an effective single impurity problem has the form of Eq. (jSJ) 
on each step of DMFT procedure. Secondly, during each DMFT iteration "external" 
self-energy Ek(zo;) is recalculated within some (approximate) scheme, taking into account 
interaction with collective modes (phonons, magnons etc.) or with fluctuations of some 
order parameter. To define non-local contribution £k(zu;) it is convenient to introduce 
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which plays the role of "bare" Green's function to build perturbation theory over "exter- 
nal" interaction. The choice of the "bare" Green's function in the form of Eq. (TTTJ) guaran- 
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"skeleton" diagrams for H^(iuj), coincides exactly with full Green's function G^(iui). 

Remarkable feature of our approach is the possibility of its generalization to calculate 
two-particle properties e.g. optical conductivity [26j [27]. Conductivity of a system is 
expressed via retarded density-density response function x R ( C0 : q) 



ie 2 uj 



-lim^-X (w,q) 



(12) 



where e is electron charge. This response function is defined by analytical continuation to 
real frequencies of the full polarization loop in Matsubara representation [28]. Note that 
conductivity is completely defined by first derivative of this response function with respect 
to q 2 in the limit of q — > 0. This circumstance, as well as the neglect of interference between 
Hubbard and "external" interactions in DMFT+S approach and locality of irreducible 
vertices of Hubbard interaction allows one to perform a partial resummation of diagrams 
relevant for conductivity, making the use of an exact (in the limit of q — > 0) Ward identity. 
At the end the real part of optical conductivity in the DMFT+S approach is [261 EZ]: 



Rea(u) 



e 2 u ro ° 



2tt 



&[/(5_)-/( £+ )]ReUf A M 



1 - 



R( 



(13) 



where 



u) = lim 



^ORR(RA) ( 



w,q) -$ 



ORR(RA) 



u,0) 



(14) 



and we introduced the two-particle Green functions of the following form: 

$° e RR{RA \uj,q) = ^G i? (e + ,k + )G^(£_,k_)r^ A )(£_,k_;e + ,k + ), (15) 



which are diagrammatically represented by Fig. [5] ( k± = k ± |, e± = e ± ^). Vertices 
T rr ( ra >(e-, k_; e + , k + ) contain all vertex corrections from "external" interaction (order 
parameter fluctuations, impurities, phonons etc.) but do not contain vertex corrections 
from Hubbard interaction. 
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Figure 5: Diagrammatic representation of <&® ra (oj, q) and $^ R (w,q). 

Thus one achieves a significant simplification of the problem. To calculate optical 
conductivity within DMFT+S approach we have only to solve single-particle problem 
of obtaining the local self-energy £(e±) with the help of DMFT+S procedure described 
above, while non-trivial contribution from non-local correlations or "external" pertur- 
bations enters via "blocks" ( Fl4|) . which can be calculated in any suitable approximation 
accounting only for "external" interaction, with "bare" Green's functions (fill , which con- 
tains the local self-energy from DMFT+S procedure. In fact Eq. ({TBI provides also an 
effective algorithm to compute optical conductivity in framework of conventional DMFT 
(neglecting all "external" non-local correlations). In this case ( I14p could be easily found 
from the simple loop diagram, defined by two Green's functions and free scalar vertices. 
To get optical conductivity there is actually no need to calculate vertex corrections in the 
framework of DMFT itself, as it was first shown for the loop with vector vertices in Refs. 

In the following, in Sections [3] and H] we shall discuss in details some applications of 
the generalized DMFT+S approach to the solution of concrete physical problems. 

2.3 Some other generalizations of DMFT 

Up to now, a number of different theoretical approaches were suggested to account for 
non-local effects within generalizations of DMFT. Here we restrict ourselves to a brief 
review of some of these approaches. 

First of all we shall refer to cluster methods already mentioned above. Instead of an iso- 
lated Anderson impurity one can consider some cluster containing several such impurities, 
treating single site correlations within DMFT, while considering intersite correlations by 
some other method. This is basic idea of the so called cluster DMFT methods [H |30| [31] . 

A drawback of such methods is related to different treatment of non-local correlation 
inside cluster and between clusters, though from the physical point of view (for example 
because of translational invariance) they should be identical. To overcome this difficulty it 
was proposed to average self-consistent "cluster" self-energy over pairs of sites connected 
by translation vectors [30]. Unfortunately such averaging procedure does not work well 
within the sel-consistent cycle of DMFT, since it breaks down analytical properties of 
Green's function. 

In some sense alternative approach, named dynamical cluster approximation (DCA), 
was proposed in Refs. [31] and J32j [33] . This approach conserves translational invariance 
and provides physical behavior of Green's function. Within DCA the Brillouin zone is 
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divided into several cells with centers denned by appropriate vectors K. The same time 
self-energy is assumed to be constant in each cell £ k (u;) = S K (w), but these constant are 
different for different K. Distinction between DCA and cluster DMFT is that in DCA 
DMFT-cluster satisfies periodic boundary conditions, instead of open boundary conditions 
in cluster DMFT. 

The choice of particular cluster method is dictated by physical problem under con- 
sideration. Naturally, the common difficulty of cluster approaches is essentially larger 
computer time consumption in contrast to the standard DMFT, which is connected with 
the solution of the appropriate cluster problem. Nevertheless, a number of successes were 
achieved on this way. Cluster DMFT generalizations were applied to different models as 
well as to studies of some realistic systems. Detailed review of these works can be found 
in Ref. [21]. 

Recently a number of diagrammatic DMFT generalizations was proposed, attempting 
for more or less consistent construction of perturbation theory over the inverse powers of 
spatial dimensionality, with the standard DMFT used as the zeroth order approximation. 
We mention in this connection the Ref. [34] and also the so called dynamical vertex 
approximation (DrA), developed in Ref. [35] . Most promising in this respect seems to 
the dual fermion approach formulated in Refs. [361137], which is claimed to be a consistent 
realization of such perturbation theory. Unfortunately up to now only some simple model 
problems were solved by these methods, while realistic systems where not yet considered 
at all. 

3 Application of generalized DMFT+E approach to 
model problems. 

3.1 Strongly correlated systems in the pseudogap state. 
3.1.1 Pseudogap fluctuations. 

Striking example of strongly correlated systems are high-T c cuprates. Parent stoichio- 
metric cuprate compounds are antiferromagnetic insulators with well developed optical 
gap and antiferromagnetism due to spin ordering on copper ions with Neel temperature 
of the order of hundreds of K. This insulating state is rapidly destroyed by introduction 
of rather few doping impurities. Thus these systems can be classified as doped Mott 
insulators with strong electronic correlations. 

Among many anomalies of the normal phase of high temperature superconductors spe- 
cial interest attracts observations of a pseudogap in the electronic spectra of underdoped 
cuprates [38, 39] . Despite continuing discussions about pseudogap nature, from our point 
of view, most preferable is the scenario of pseudogap formation due strong scattering of 
charge carriers on antiferromagnetic (AFM, SDW) short range order fluctuations [39] 140]. 
In momentum representation this scattering is characterized by momentum transfer of the 
order of Q = (-, -) (a is two-dimensional lattice constant). This leads to formation of 
certain features in single-particle spectrum, which are precursors of changes in the spectra 
due to the appearance of AFM long range order (period doubling). As a result we end 
up with non Fermi liquid behavior (dielectrization) of spectral density in the vicinity of 
the so called "hot-spots" on the Fermi surface, appearing at intersections of Fermi surface 
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with borders of AFM Brillouin zone 

In the framework of this spin-fluctuation scenario in works a simplified model of the 
pseudogap state was studied in Refs. [391 El H2] ■ This model is based on the assumption 
that for high enough temperatures dynamics of spin fluctuations can be neglected and 
one can consider instead the scattering of charge carriers by static Gaussian random 
field (quenched disorder) of pseudogap fluctuations (short range order AFM fluctuations). 
Scattering intensity on fluctuations is characterized by a narrow peak near scattering 
vectors of the order of Q with a width defined by inverse correlation length k = £ _1 and 
corresponding energy scale A (of the order of pseudogap crossover temperature T*). 

Thus for momentum dependent self-energy we shall concentrate on the case of electron 
scattering on such (SDW-like) antiferromagnetic spin fluctuations (similar consideration 
works fine also for CDW-like charge fluctuations) with short range order. To calculate 
Sk(icu) for the case of electrons propagating in quenched random field of Gaussian spin (or 
charge) fluctuations with dominating scattering momentum close to characteristic vector 
Q ("hot-spot" model [2H]), we shall use the generalized version of recurrent procedure 
proposed in Refs. [HI H21 H3] , allowing to take into account all Feynman diagrams de- 
scribing scattering of electrons by this random field. This becomes possible because of 
remarkable property of simplified "hot-spot" model: contribution of arbitrary diagram 
with crossing interaction lines is equal to contribution of some diagram of the same order 
without crossing of those lines [13] . Thereby we can restrict ourselves to consideration of 
non-crossing diagrams only and take into account other diagrams by combinatorial pref- 
actors attributed to interaction lines [121 H3] . Finally we obtain the following recurrent 
relation for the self-energy (continuous fraction representation [32l 133] ): 

EJiuk) = A 2 — — . (16) 

icj + fi — Li[ico) — £ n (kJ + mv n K — h n+ i{iu, k) 

Here, the term E n (iu;,k) of recurrent procedure contains all diagrammatic contributions 
with number of interaction lines > n. Recurrent procedure for T, n (iu, k) converges rather 
fast, we can put S n (icu,k) for large enough n equal to zero and performing straightfor- 
ward computations obtain the desired physical self-energy for n = 1 [35] , which can be 
subsequently used in DMFT+E computational scheme: 

S k (iw) = £ n=1 (iu;,k) (17) 

Parameter A characterizes the energy scale (width) of the pseudogap, k = is 
the inverse correlation length of SDW (CDW) fluctuations, e n (k) = e(k + Q) and v n = 
l^k+cjl + l^k+ol f° r °dd n, e n (k) = e(k) and v n = \v£\ + \v^\ for even n, where velocities 
projections t>£ and are defined by usual derivatives with respect to corresponding 
momenta components of the bare electron dispersion e(k). At last, s(n) are combinatorial 
prefactors defining the number of coinciding diagrams: 

s(n) = n (18) 

for the case of commensurate charge (CDW-type) fluctuations with Q = (n/a,ii/a) [13]. 
For incommensurate CDW fluctuations [33] (when Q is not related to lattice period) we 
get: 

/ \ / r± r for odd n f . 

^ ^ for even n. 
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If we want to take into account spin (Heisenberg) structure of interaction with spin fluc- 
tuations for nearly antiferromagnetic Fermi liquid (spin-fermion model [H|), the combi- 
natorics becomes more complicated. Scattering processes preserving spin projection are 
controlled by commensurate combinatorics, while spin flip scattering is described by dia- 
grams of incommensurate type ("charged" random field, according to Ref. [Hj). In this 
model the recurrent procedure (fT6j) for single-particle Green's function remains the same 
but with another combinatorial pref actors s(n) |JT]: 

S (n) = { fOT ° dd H (20) 

\ j for even n. 

Obviously this procedure introduces an important length scale £, missed in standard 
DMFT. Physically this length scale reflects the influence of short range order fluctuations 
(SDW or CDW) on electronic bath surrounding effective Anderson impurity in DMFT. 

After the self-consistent solution of DMFT+S set of equations (l5l fT0l) is obtained, one 
can calculate the spectral density A(oj,k): 

A(u,k) = —1m 7Z r 1 w r = . r , (21) 

v ' 7i u + n - e(k) - E(w) - £ k (u;) ' v ' 

where £(u;), Xk(w) and chemical potential fi are already computed in a self-consistent 
way. Density of states can be found by integration of (12 II) over the Brillouin zone. 

Analogous approach can be developed also to determine two-particle vertices. Basic 
idea employed here is the possibility to get arbitrary vertex diagram by introducing "ex- 
ternal field" line into corresponding self-energy diagram [44j |45j 146]. In the model under 
consideration we can again restrict ourselves to non-crossing diagrams, while contribution 
of all other diagrams can be accounted for by combinatorial prefactors s(n) attributed 
to interaction lines [^T| H2J H3] . Thus, all vertex diagrams are obtained from the simple 
ladder diagrams with additional prefactors s(n) on corresponding interaction lines [15|H6] 
(see also [29J). Then we obtain the following system of recurrent relations for the vertex 
T ra (e_, k_; e + , k + ) [46J, where contribution of local DMFT self-energy (obtained within 
the DMFT+E procedure) is already included: 



rj* (e_,k_;e + ,k + ) = 1 + A 2 S (n)G£(e_, k_)G£ (e +l k^ 



x 



2,iv n Kk 

u - e n (k + ) + e n (k_) - + S^(e_) - Z* +1 (e + , k + ) + S^ +1 ( £ _, k_) 

xr^(e_,k_;e + ,k + ), (22) 

where 

G*> A (e±, k±) = g-i . (23) 

£±-£ n (k ± )±Z W -E^( £± )-E^(£ ± ,k ± ) 

"Physical" vertex T RA (e-, k_; e + , k + ) is defined as r^f (£_, k_; e + , k + ). Recurrent pro- 
cedure ( 1221) accounts for a// diagrams of perturbation theory for the vertex part. In the 
limit of k — y (£ — >• oo) f l22|) can be reduced to a series investigated in Ref. [S] (see also 
[H]), which can be exactly summed in analytical form. Standard ladder approximation 
is reproduced if all combinatorial factors in ( 122]) are made equal to one for all n |45j . 
Recurrent procedure for T RR (e^, k + ; e + , k + ) differs from (j22l only by the evident change 
of A — > i?, as well as replacing the whole expression in figure brackets on the r.h.s. of 
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Eq. (J22D by 1. Eqs. OH), ([IS]), 022]) together with flU) and OH provide the complete self- 
consistent procedure to calculate optical conductivity within our model in the framework 
of DMFT+S approach. 

Important aspect of our theory is the possibility of microscopic calculation of both 
effective parameters A and £. For example, applying two-particle self-consistent theory of 
Ref. |47J, together with approximations introduced in Refs. [2J H2] for two-dimensional 
Hubbard model, we derived a microscopic expression for A [23], which can be calculated 
within the standard DMFT. It can be shown that for wide range of hole doping the 
pseudogap amplitude A varies in the interval from t to 2t (t is the nearest neighbor 
hopping integral). 

3.1.2 Basic electronic properties in the pseudogap state. 

Let us discuss results for the standard single band Hubbard model on a square lattice 
with electron dispersion 

e(k) = — 2i(cos k x a + cos k y a) — At' cos k x a cos k y a , (24) 

with t and t' nearest and next nearest hopping integrals. 

Energy scale in the following is defined by nearest neighbor hopping integral t, and 
length scale by the lattice constant a. Impurity solver used was the numerical renormal- 
ization group (NRG) [HI EES]- Detailed computational results on single particle properties 
demonstrating pseudogap anomalies can be found in Refs. [231 [2H 125], and on optical 
conductivity in Ref. [26] . Here we only discuss most typical results corresponding mostly 
to the case of t'/t = —0.4 (characteristic for cuprates) and band filling n = 0.8 (hole 
doping) . 

Density of states and spectral function. 

Lets start with results obtained within generalized DMFT+S approach for the densi- 
ties of states (DOS) in case of rather weak (compared to bandwidth) Coulomb interaction 
U = At. Characteristic feature of strongly correlated metallic state is coexistence of lower 
and upper Hubbard subbands splitted by Coulomb interaction U with quasiparticle peak 
at the Fermi level [U [5J. Noninteracting DOS for the square lattice has Van- Hove sin- 
gularity near the Fermi level, so that in general the peak on the Fermi level can not be 
treated simply as a quasiparticle one. Actually there are two contributions to this peak: 
(i) from quasiparticle peak appearing in strongly correlated metals because of manybody 
effects and (ii) smoothed Van- Hove singularity of noninteracting DOS q 

On the left side of Fig. [6] we show DMFT+S DOS with n = 0.8 for the case of 
t'/t = —0.4 (left panel) and t' = (right panel) for two different temperatures T = 
0.356t (lower panel) and T = 0.088t (upper panel). Black curves are obtained in the 
absence of fluctuations. Other curves on Fig. [6] present results for DOS with nonlocal 
fluctuations with amplitude A = It. For all parameter sets we see that introduction 
of non-local fluctuations leads to pseudogap formation on a quasiparticle peak. Width 
of the pseudogap (energy interval between corresponding peaks in the DOS) is of the 
order of ~ 2 A. Decrease of the value of A from 2t to t gives twice smaller pseudogap 
width and makes it less deep. More pronounced pseudogap is found for spin-fermion 

3 With decrease of Coulomb repulsion Van-Hove singularity gradually transforms into quasiparticle 
peak at U = (6 + 8)i. 
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Figure 6: Comparison of DMFT(NRG)+£ DOS calculated in Ref. [23] for different combi- 
natorial prefactors (SF — spin-fermion model, commensurate fluctuations), inverse correlation 
lengths £ _1 (in units of lattice constant), pseudogap amplitude A = 2t and band filling n = 0.8. 
On the left side — U = 4t, t'/t = —0.4 (left column), t' = (right column), temperature 
T = 0.088t (upper part) and T = 0.356t (lower part). On the right side -U = 40t, T = 0.088*. 
Fermi level corresponds to zero energy. 

model combinatorial prefactors (see fl20l) ). as compared with the case of commensurate 
charge fluctuations (combinatorial prefactors (|T9|) ). The influence of correlation length 
value corresponds to general expectations. Decrease of correlation length or, for inverse 
correlation length, the change from = 0.1 to = 0.5 slightly smears the pseudogap. 
The rise of temperature from T = 0.088t to T = 0.356t leads to general broadening of 
DOS structures. One should note that DMFT+S results for U = 4t (which is less than 
the bandwidth W) are qualitatively similar to results obtained earlier in the absence of 
Hubbard interaction [HI [42] . 

Let us consider now the case of doped Mott insulator with Hubbard interaction value 
U = 40t, t'/t = —0.4 and band filling n = 0.8 (right side of Fig. [6]). Characteristic feature 
of DOS for such strongly correlated metals is strong splitting of lower and upper Hubbard 
bands with the Fermi level within the lower Hubbard band (the case of hole doping) In 
absence of nonlocal fluctuations again the quasiparticle peak is formed on the Fermi level. 
However upper Hubbard band now is quite far away to the right and does not touch the 
quasiparticle peak (as it does for the case of weak Hubbard interaction). 

For strong enough nonlocal fluctuations with A = 2t pseudogap appears in the middle 
of quasiparticle peak and the lower Hubbard band is slightly broadened by fluctuations 
effects. Qualitatively pseudogap anomalies behavior reminds that described above for 
the case of U — At - decrease of £ smears the pseudogap and makes it less pronounced, 
decrease of A from A = 2t to A = t narrows the pseudogap and makes it more shallow 
(see. |23j). Let us notice also that for the doped Mott insulator pseudogap is more evident 
for spin SDW-like fluctuations than for the charge CDW-like ones. 

Nevertheless there are quite appreciable distinctions in contrast to the U = 4t case. 
For example, the width of the pseudogap in DOS is found to be essentially smaller than 
2A which is connected, in our opinion, with noticeable narrowing of the quasiparticle 
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Figure 7: DMFT(NRG)+S spectral densities A(k, co) [23] along high-symmetry directions of 
the first Brillouin zone r(0, 0) — X(7r, 0) — M(7r, tt) — T(0, 0), for spin-fermion combinatorics (SF). 
The Fermi level is at zero energy. 

peak itself caused by local correlations. 

In Fig. [7] we show spectral densities A(u,k), calculated within the DMFT+E ap- 
proach along high-symmetry directions of the first Brillouin zone: T(0, 0) — X(7r, 0) — 
M(7r, 7r) — r(0, 0). In fact this figure shows the quasiparticle band of manybody system 
- positions of maxima of spectral functions define quasiparticle dispersion, while their 
width defines quasiparticle damping. Also we clearly observe the partial reconstruction 
( "destruction" ) of this band by pseudogap fluctuations. One can see characteristic double- 
peak structure close to X-point of the Brillouin zone. In the middle of M—T direction (the 
so called "nodal"-point) one can observe the rise of the pseudogap, i.e. the "memory" of 
the AFM gap, which has maximum here in the case of AFM long range order. Generally 
speaking varying of filling leads to a shift of spectral functions with respect to the Fermi 
level. 

Fermi surface "destruction". 

Within conventional DMFT Fermi surface is not renormalized by interaction i.e. it 
stays the same as for the bare quasiparticles [3J. However in the case of nontrivial self- 
energy momentum dependence substantial renormalization of the Fermi surface appears 
due to pseudogap formation [JT]. There are several ways to define Fermi surface for 
strongly correlated systems. Below we shall exploit intensity map of the spectral function 
(spectral density) ( 12 ip for u — 0, which is often called the Fermi surface map. Such a map 
is directly measured by ARPES experiments and positions of its intensity maxima specify 
the Fermi surface in a sense of the usual Fermi liquid theory, in case of quasiparticle 
damping becoming negligibly small. 

On the right side of the Fig. [8] there are displayed these maps for case of strongly cor- 
related metal with U = At. This contour plot of the spectral function clearly demonstrates 
the "destruction" of the Fermi surface in the "hot-spots" together with formation of the 
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Figure 8: DMFT(NRG)+S picture of the Fermi surface "destruction" Ref. [Mj On the left - 
for U = At. On the right - for U = 40t (a) A = 0.2*; (b) A = 0.4t; (c) A = t; (d) A = 2t. Band 
filling is n = 0.8. Dashed line - "bare" Fermi surface. Full lines - solutions of Eq. (125 ^ . 

"Fermi arcs" upon A growth similar to that observed in pioneer works of Norman et al. 
[4"8] , which were later confirmed in a large number of other works. One should note that 
qualitatively analogous behavior is also found in the absence of local electronic correlations 
(U = 0) pU [49]. The role of finite U values add up to the decrease of spectral function 
intensity compared to the case U = and leads to additional "smearing" making "hot- 
spots" less visible. "Destruction" of the Fermi surface starts in the vicinity of "hot-spots" 
for small A values, but practically simultaneously it disappears in the whole antinodal 
region (near points X(7r,0), Y(0,7t)) of the Brillouin zone, while only "Fermi arcs" in the 
nodal region remain, with the shape close to bare Fermi surface. Those results naturally 
explain why in ARPES the clear "hot-spots" behavior is rather rarely observed [50]. In 
more details the question of possibility to observe "hot-spots" will be elucidated in the 
section below, devoted to LDA+DMFT+S description of realistic cuprates. 

In the case of doped Mott insulator with U = 40t shown in Fig. [8] we see that "Fermi 
surface" is rather ill defined for all values of A. The profile of spectral function is sig- 
nificantly more "smeared" in contrast to smaller values of U reflecting important role of 
local correlations. For comparison in Fig. E]we also show the renormalized Fermi surfaces 
obtained within the model by formal solution of Eq. ( 125]) 0: 

w-e(k) +fi-ReZ(cu) - i?eS k (w) = (25) 

for uj = used for example in the work [H]. Obviously, this definition gives Fermi 
surface close to the one obtained from intensity map for small A, but does not account 
for significant damping essentially important for large A. For large pseudogap amplitudes 
this definition of the Fermi surface is qualitatively adequate to true behavior, resulting 

4 This definition is used in the standard Fermi liquid theory. In fact in our particular case influence of 
nonlocal pseudogap fluctuations leads to qualitative changes of simple Fermi liquid picture. Herewith we 
leave aside the question whether Fermi liquid is applicable for such defined Fermi surface or not in the 
limit T — > since static approximation used is by construction the high temperature one - short range 
order AFM fluctuations can be considered as quasistatic only if T ^> u s f, where tu s f - characteristic 
frequency of spin fluctuations j41] [42] . 
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Figure 9: Real part of DMFT+E optical conductivity (t! = -0.4t, t = 0.25 eV) obtained in 
Ref. [26] for different values of the pseudogap amplitude: A = 0, A = t, A = 2t . Temperature 
is T = 0.088t, band filling - n = 0.8 and correlation length £ = 10a. On the left side - strongly 
correlated metal with U = At. On the right side - doped Mott insulator with U = 40i. Inset: 
conductivity in a wide frequency range, which includes transitions to the upper Hubbard band. 

from spectral function analysis, only in the nodal region. Actually, the contour plot of 
spectral function (at u — 0) gives most complete and natural representation of the Fermi 
surface for the systems with strong correlations and nonlocal fluctuations of some order 
parameter, which are present in a wide region of the phase diagram of high-T c cuprates 
because of their low dimensionality. Results obtained in a such an approach directly 
correspond to ARPES experiments, where exactly this definition of the Fermi surface is 
most conventional. 

Optical conductivity. 

Lets set about discussion of DMFT+S results for optical conductivity. On the left 
panel of Fig. [9] we show DMFT+ S results for the real part of optical conductivity in the 
case of strongly correlated metal (U = At) for different values of the pseudogap amplitude. 
We clearly observe the formation of typical pseudogap anomaly on the "shoulder" of the 
Drude peak and it grows as A increases. This behavior is rather similar to "mid-infrared 
feature" which is observed in optical conductivity of cuprate superconductors j5TJ[52]. The 
rise of temperature and decrease of fluctuations correlation length wash off pseudogap, 
making this anomaly less pronounced [26J. 

The right panel of Fig. [9] demonstrates DMFT+ S optical conductivity of doped 
Mott insulator (U = AOt) for several values of the pseudogap amplitude. We see that 
frequency range where pseudogap anomaly is observed gets narrower with the growth 
of local correlation strength and for large U values pseudogap anomalies are strongly 
suppressed. Pseudogap fluctuations lead to noticeable changes of optical conductivity 
only for relatively low frequencies, of the order of A. For higher frequencies (e.g. of the 
order of U, where transitions to upper Hubbard band take place) pseudogap effects are 
not seen (see also inset on the right panel of Fig. [9]). For low frequencies we observe 
suppression of Drude peak with rather weak anomaly at u> ~ A which disappears for 
small A values or for short correlation lengths. 
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3.2 Mott- Anderson transition in disordered systems. 

The importance of both electron interactions and disorder effects in the condensed matter 
research is well known [53J. Coulomb interaction and disorder are two driving forces 
leading to metal-insulator transition, connected with localization and derealization of 
charge carriers. In particular Mott-Hubbard transition is induced by electron repulsion 
[H EI] , while Anderson metal-insulator transition is related to scattering of noninteracting 
particles by impurities [55]. It is well known that a subtle competition between disorder 
effects and interaction has many manifestations [53~1 [56] . Most relevant this problem 
is in the case of strong disorder and strong electron correlations, determining physical 
mechanisms of Mott- Anderson metal- insulator transition [53] . 

One of the main models allowing for the account of both electronic correlations (leading 
to Mott metal- insulator transition [UEI]) an d strong disorder effects (leading to Anderson 
metal- insulator transition) is Anderson-Hubbard model [27J E7J ESI Ell EH ED E2] • 

In Refs. J57J ESI EH] three-dimensional Hubbard-Anderson model was investigated 
in the framework of dynamical mean-field theory (DMFT) [3 El El II]- Influence of lo- 
cal disorder was taken into account through averaged density of states (DOS) [631 El] 
within the well known coherent potential approximation (CPA), which does not describe 
Anderson localization. To overcome this difficulty in Dobrosavljevic and Kotliar [57] 
has proposed a version of DMFT, where the self-consistent solution of stochastic DMFT 
equations for an ensemble of systems with given realizations of disorder, was used to 
calculate the averaged logarithmic (geometric mean) density of states, which gives infor- 
mation on critical disorder for Anderson transition. Further this approach was developed 
in Refs. [58] [59] where highly nontrivial phase diagram of three-dimensional paramagnetic 
Anderson-Hubbard model [59] was obtained, containing correlated metal phase, Mott in- 
sulator phase and correlated Anderson insulator phase. Main problem of the approach 
used in Refs. [571 ESI EH] is the impossibility of direct computation of measurable physical 
properties such as conductivity, which actually defines metal-insulator transition. 

At the same time there exists the well developed self-consistent theory of Anderson 
localization, based on the solution of equations for the generalized diffusion coefficient. 
Efficiency of this approach in the absence of interactions is known for a long time [2S1 ESI 
I3H] IB"5] IB"B"1 157] , certain attempts to include interaction effects into this approach with some 
promising results were undertaken in Refs. [661 ES] ■ However, up to now this approach 
was not extended to modern theory of strongly correlated systems. For the first time 
such investigation was performed in Ref. [27] for three-dimensional systems and later for 
two-dimensional case [69"] . 

Let us consider disordered paramagnetic Anderson-Hubbard model (mostly) at half- 
filling for arbitrary interaction and disorder strength. Obviously this model contains both 
Mott-Hubbard and Anderson metal-insulator transitions. Hamiltonian of the model is: 



where t > is nearest neighbor hopping amplitude, while U is on-site Hubbard repulsion, 
n ia = a la a ia * s particle number operator, a ia (a\ a ) is annihilation (creation) operator of 
electron on site i with spin a. Local energies e« are assumed to be random and independent 
at different lattice sites. To simplify diagram technique hereafter we assume the Gaussian 




(26) 
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distribution for e 



y^A V 2A 2 

Parameter A here is the measure of disorder. Such Gaussian random field ( "white noise" ) 
of energy levels e» at different lattice sites is equivalent to impurity scattering and can be 
described by standard diagram technique for the averaged Green functions [29J. 

Self-energy E p (ze) caused by scattering by disorder calculated in simple one- loop ap- 
proximation neglecting "crossing" diagrams (i.e. in self-consistent Born approximation) 
[29] , in case the of Gaussian disorder (1271) reduces to: 



E p (z £ ) = A 2 J2 G{ie, p) = S imp (ze), (2* 



so that our "external" self-energy is independent of momentum p (local). 

To analyze optical conductivity we shall apply the general DMFT+S expression ([TBI . 
Most important block $° RA (oj,q) can be obtained using the ideology of self-consistent 
theory localization [281 12HI ESI [39j [661 EI], with some generalizations to account for the 
role of Hubbard interaction via DMFT+S approach [27J EH]. Main distinction from 
the standard derivation of the equations of self-consistent theory of localization is the 
use of Green's functions (j4j) containing local contribution to self-energy from Hubbard 
interaction. 

Following standard derivation [281 [29J [391 ESI EHl EZ] we obtain diffusion-like (for small 
oj and q) contribution to $° (u), q) which takes the form: 

u + iD(oj)q z 

where D(u) is the generalized diffusion coefficient and important difference from the 
single-particle case is contained in: 

Q = s + -s_ - Z R (s + ) + S A (e_) =w- £*(£+) + Z A (s_) =w- A£ RA (u;), (30) 

which substitutes for the usual u term in the denominator of standard expression for 
$® ra (cj, q). From general considerations it is clear that in metallic phase for uj — > we 
have AE M (w = 0) = 2ilmT,(e) ~ Max{T 2 ,e 2 }, which reflects Fermi liquid behavior in 
DMFT (which is not violated by elastic impurity scattering). For finite T it leads to usual 
phase decoherence caused by (inelastic) electron-electron scattering [531 156]- 
Then Eq. f[T3|) takes the form: 

AE rr (oj 



. e 2 co f°° , rr . . P . ) 27TN(e)D(uj) , nRR/ , 
Rea(u) = — J de [/(e_) - /(e + )] Re j - $ rr (uj) 

(31) 

where the second term in figure brackets can be in fact neglected for small u, while in 
case of necessity to describe a wide frequency range it can be calculated using ( TH| with 
§® RR (u, q) taken in the usual ladder approximation. 

Now following the standard scheme of self-consistent theory of localization [2H1 EHl EH 
ESI E51 E7J we get the closed self-consistent equation for generalized diffusion coefficient: 

DH = - AE- H + A^(A Gp )^ -^-^J (32, 
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where d is spatial dimensionality, AG P = G R (e + , p) — G A (e_, p), A£^£(u;) = p (s+) — 
(eJ) and averaged velocity < v >, which can be well approximated just by the Fermi 
velocity, is given by the following expression: 

V IvJAGp 

< ' >= IzzgT' (33) 

where v p = Eq. (!32|) should be solved together with self-consistent DMFT+E 

procedure floTfTUl |2"B"]) . In fact, this equation is a transcendental one and can be easily 
solved by iterations for each Q value. 

In accordance with the usual applicability limit of diffusion approximation, summation 
over q in Eq. (132]) should be restricted to [66, [29] : 

g<fc = Mm{r 1 ,p F }, (34) 

where / =< v > /2 r y is an elastic mean free path, 7 is Born scattering frequency by 
impurities, Pf is the Fermi momentum. It is well known that in two-dimensional case 
Anderson localization occurs for any infinitely weak disorder. However, in this case the 
localization radius is exponentially large and sample size becomes essentially important. 
Sample size L can be introduced into self-consistent localization theory introducing the 
integration cutoff of diffusion pole at small q [28, 65J, i.e. at: 



q ~ k L = l/L. (35) 

For oj (on the Fermi surface (e = 0), and obviously Co — > 0) in the Anderson 
insulator phase one gets localization behavior of the generalized diffusion coefficient [28] 
65] [29]: 

D{u) = -luRioc- (36) 
After substitution of f[36]) into f[32l one obtains equation defining localization radius Ri oc '- 



Rloc 



< v > 
rfA 4 



2 



3.2.1 Three-dimensional systems. 

Below we present most interesting results for three-dimensional Anderson-Hubbard model 
at half-filling on a cubic lattice with semielliptic bare density of states with the bandwidth 
W = 2D: 

N °( £ ) = ^i^^- ( 38 ) 

Density of states is given in units of number of states in energy interval for unit cell of the 
volume a 3 (a is the lattice constant) and for one spin projection. Conductivity values are 
given in natural units of e 2 /ha (a is the lattice constant). For more detailed acquaintance 
with numerical results (also for the case of deviations from half-filling) we refer the reader 
toRef. [27]. 

Evolution of the density of states 

Within the standard DMFT approach Hubbard model density of states at half-filling 
has a typical three-peak structure [5] H] [70] with a narrow quasiparticle peak (central 
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Figure 10: Hubbard- Anderson model density of states at half- filling for different disorder levels 
A [27]. On the left side — correlated metal with U = 2.5D, On the right side — Mott insulator 
with U = 4.5D. 

peak) at the Fermi level and wide upper and lower Hubbard bands situated at energies e ~ 
±U /2. AsU grows quasiparticle band narrows within the metallic phase and disappears at 
Mott-Hubbard metal-insulator transition at critical interaction value U C 2 ~ 1.5VF. With 
further increase of U insulating gap opens at the Fermi level. 

In Fig. [10] we present our results for DMFT+E densities of states for typical strongly 
correlated metal with U = 2.5D = 1.25W, in the absence of disorder and for different 
values of disorder A, including strong enough disorder, transforming correlated metal to 
correlated Anderson insulator (see also next section on conductivity). As one can expect 
disorder leads to typical broadening and suppression of the density of states. 

More unexpected is the result obtained for U = 4.5D = 2.25W, typical for Mott insu- 
lator and shown on right panel of Fig. [TJJ Here we observe restoration of the central peak 
(quasiparticle band) in DOS with the increase of disorder, transforming Mott insulator to 
correlated metal or to correlated Anderson insulator. Similar DOS behavior was reported 
also in Ref. [59). 

Physical origin of such quite unexpected central peak restoration is pretty clear. Con- 
trolling parameter of metal-insulator transition in DMFT is the ratio of Hubbard inter- 
action U to bare bandwidth W = 2D. With disordering (in the absence of Hubbard in- 
teraction) new effective bandwidth W e ff appears which grows with disorder. Semielliptic 
form of DOS with well defined band edges within the self-consistent Born approximation 
( l28l) is preserved. This leads to diminishing values of the ratio U/W e ff, which in its turn 
causes restoration of the quasiparticle band. In more details this is discussed below, when 
we discuss the phase diagram of Hubbard-Anderson model. 

In absence of disorder characteristic feature of Mott-Hubbard metal-insulator transi- 
tion is hysteresis DOS behavior appearing with the decrease of U starting from insulating 
phase [5] [70]. Mott insulator phase is conserved (metastable) down to rather small U 
values deep within the correlated metal phase. Metallic phase is restored only at about 
U c \ ~ LOW. Corresponding interval U& < U < U& typically is considered as a coexis- 
tence region of metallic and Mott insulating phases, where, from a thermodynamic point 
of view, metallic phase is more stable [5] [701 [71] • Such hysteresis DOS behavior (see 
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co/2D 

Figure 11: Real part of optical conductivity of Hubbard- Anderson model at half- filling for 
different disorder levels A [27J. On the left side - typical correlated metal with U = 2.5D. Curves 
1,2, - metallic phase, curve 3 corresponds to the mobility edge (Anderson transition), curves 
4,5 - correlated Anderson insulator. On the right side typical Mott insulator with U = 4.5D. 
Curves 1,2 correspond to Mott insulator, curve 3 - mobility edge (Anderson transition), curves 
4,5 - correlated Anderson insulator. Inset - magnified low frequency region. 

[27] ) is observed also in the presence of disorder and will be described below during the 
discussion of the phase diagram of Anderson-Hubbard model. 

Optical conductivity: Mott-Hubbard and Anderson transitions 

Without disorder our calculations reproduce conventional DMFT results [H [5] , where 
optical conductivity is characterized by the usual Drude peak at low frequencies and wide 
maximum at about ou ~ U, which corresponds to optical transitions to upper Hubbard 
band. As U grows Drude peak diminishes and disappears at Mott transition. Introduction 
of disorder leads to qualitative change of the frequency dependence of optical conductivity. 

On the left panel of Fig. [11] we show the real part of optical conductivity of Hubbard- 
Anderson model at half-filling for different disorder levels A and U = 2.5D typical for 
correlated metal. Transitions to the upper Hubbard bands at energies co ~ U are almost 
unobservable. However it is clearly visible that metallic Drude peak typically centered 
at zero frequency is broadened and suppressed by disorder, gradually transforming into 
a peak at finite frequency because of Anderson localization effects. Anderson transition 
takes place at A c « 0.74D = 0.37H 7 (corresponding to the curve 3 on all figures here and 
also for DOS). Notice that this value depends on the cutoff (l34"l) . which is defined up to the 
coefficient of the order of unity [66], EH]. Naive expectations can bring us to a conclusion 
that narrow quasiparticle band at the Fermi level (formed in a strongly correlated metal) 
may be localized much easily than normal conduction band. However we see that these 
expectations are wrong and the band localizes only at rather large disorder A c ~ D, 
similar to that for conduction band of the width ~ W. It agrees with the known analysis 
of localization in a two-band model [72] . 

In the DMFT+E approach critical disorder value A c does not depend on U as inter- 
action effects enter Eq. (1321) only through AS M (u) — > for oj — > (for T = 0, e = 0), 
and the influence of interaction at u> = disappears. In fact this is the main shortcoming 
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of DMFT+S approach originating from the neglect of interference effect between inter- 
action and impurity scattering. Significant role of these interference effects is known for 
a long time [531 [56]. O n the other hand, the neglect of these effects allows to perform the 
reasonable physical interpolation between two main limits - that of Anderson transition 
because of disorder and Mott-Hubbard transition because of strong correlations. One 
can consider this approximation as a reasonable first step to a future complete theory of 
metal-insulator transition in strongly correlated disordered systems. 

On the right panel of Fig. [11] we show the real part of optical conductivity of Mott- 
Hubbard insulator with U = 4.5D at different disorder levels A. In the inset we show low 
frequency data, demonstrating different types of conductivity behavior, especially close to 
Anderson transition and within the Mott insulator phase. On the main part of the figure 
contribution to conductivity from transitions to upper Hubbard band at about u ~ U is 
distinctly seen. Disorder growth results in the rise of finite conductivity for the frequencies 
inside Mott-Hubbard gap, correlating with the restoration of quasiparticle band in DOS 
within the gap as shown in the right panel of Fig. [TD1 This conductivity for A < A c 
is metallic (finite in the static limit u — 0), and for A > A c at low frequencies we get 
Rea(u) ~ u 2 , which is typical for Anderson insulator [28] [291 ES ESI EHl IB"?] . 

A bit unusual is the appearance in ~Rea(u) of a peak at finite frequencies even in 
the metallic phase. This happens because of importance of localization effects. In the 
"ladder" approximation for Q® RA (cu,q) which neglects all localization effects we obtain 
the usual Drude peak at u = [2?], while account of localization effects shifts the peak 
in Re<7(co>) to a low (finite) frequencies. As is well known [TJ, metallic state is defined by 
finite static (ui = 0) conductivity at zero temperature. 

Above we presented the data for conductivity data obtained for the case of increase of 
U from metallic to Mott insulator phase. As U decreases from Mott insulator phase we 
observe hysteresis of conductivity in coexistence region defined (in the absence of disorder) 
by inequality JJ C \ < U < U C 2- Hysteresis of conductivity is also observed in the coexistence 
region in the presence of disorder. Details of this behavior of optical conductivity can be 
found in Ref. |27J. 

Phase diagram of Anderson-Hubbard model at half-filling 

Phase diagram of Anderson-Hubbard model at half-filling was studied in Ref. [59], 
using direct DMFT calculations for the lattice with finite number of sites with random 
realizations of energies in ( 126]) and averaging over these realizations to get averaged 
DOS and geometric mean local DOS which allows one to define critical disorder for tran- 
sition into Anderson insulator phase. Below we present our results on Anderson-Hubbard 
model phase diagram obtained from DOS and optical conductivity calculations within 
the DMFT+E approach. One should emphasize that conductivity analysis is most direct 
way to distinguish metallic and insulating phases [I]. 

Calculated disorder-correlation (A, U) phase diagram at zero temperature is shown 
on the left panel of Fig. [12] Anderson transition line A c w 0.37W = 0.74D is defined as 
a disorder strength for which static conductivity becomes zero at T = 0. Mott-Hubbard 
transition can be detected from central (quasiparticle) peak disappearance in DOS or 
from optical conductivity by observation of gap closing in the insulating phase or from 
static conductivity disappearance in the metallic phase. 

We have already noticed that DMFT+S approximation gives universal (U indepen- 
dent) value of critical disorder A c because of neglect of interference between disorder 
scattering and Hubbard interaction. This leads to differences between phase diagram of 
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Figure 12: Phase diagram of paramagnetic Anderson-Hubbard model. On the left side - 
zero temperature case |27| . Continuous curves are Mott insulator phase boundaries C/ c i jC 2(A) 
obtained from analytical estimate of Eq. (|4ip . different symbols represent results for these 
boundaries obtained from calculations from DOS and optical conductivity. Line of Anderson 
transition is given by A c = 0.37. On the right side - finite temperature case. Points are obtained 
from DOS calculations. Solid black curve is linear fit (T C /2D = 0.02(1 + A/2D)) to the T c points 
where coexistence region disappears. 



Fig. [12] and the one obtained in Ref . [59J . At the same time influence of disorder scat- 
tering on Mott-Hubbard transition is highly nontrivial and qualitatively coincide with 
results of Ref. [59] . Main difference is conservation of Hubbard bands in our results even 
in the limit of high enough disorder, while in the Ref. [59] they just disappear. Moreover 
coexistence region in Fig. [12] slowly widens with disorder growth instead of vanishing at 
some "critical" point as on phase diagram of Ref. [59]. Coexistence boundary regions, 
which are denned by Mott insulator phase boundaries, obtained with decrease on increase 
of U, represented by curves U c i(A) and U C 2(A) on Fig. [121 can be obtained from the 
simple equation: 

where effective bandwidth in the presence of disorder was calculated for U = within 
self-consistent Born approximation (|28|) : 



W eff = W^ 1 + 16^. (40) 

Thus the boundaries of coexistence region which define also Mott insulator phase bound- 
aries are given by: 

U cliC2 (A) = U cliC2 ^ 1 + 16^ (41) 

which are shown in Fig. [12]by dotted and solid lines. Phase transition points detected from 
disappearance of quasiparticle peak as well as points following from qualitative changes 
of conductivity behavior are shown in Fig. [12] by different symbols. These symbols 
demonstrate very good agreement with analytical results supporting the choice of ratio 
as a control parameter of Mott transition in presence of disorder. 
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On the right panel of Fig. [T2]we show temperature dependence of Mott insulator phase 
boundaries U C 2(A) (squares) and coexistence region U cl (A) (circles) @. It is well known 
that in "pure" DMFT without disorder coexistence region narrows with temperature 
growth and vanishes at some critical temperature T c . In the DMFT+S approach disorder 
increase at zero temperature widens the coexistence region and this behavior remains 
for finite temperatures. It is also seen that disorder growth leads to a practically linear 
growth (solid black line) of this critical temperature. Notice also very weak temperature 
dependence of coexistence boundary C/ cl (A). 

3.2.2 Two-dimensional systems. 

According to scaling theory of localization [73] metallic state in two-dimensional (2D) 
systems does not exists, electrons are localized already at any infinitely small disorder. 
Despite this prediction for 2D systems was made for noninteracting particles, later it was 
found that in the simplest case weak interaction between electrons also favors localization 
[73] . In early 80th experiments done on different 2D systems [75] mostly confirmed these 
predictions. However, some theoretical works [56] pointed that this point of view in 
general is incorrect, since in the limits of weak disorder and large enough interaction 
2D systems can have finite conductivity at zero temperature. Experimental discovery of 
metal-insulator transition in 2D weakly disordered systems at low carrier concentration, 
absent in the single particle theory, stimulated new direction of theoretical studies (see 
introduction to this field in review papers [76| 177]). 

In DMFT+E approximation, as we shall see, for infinite 2D system (L — > oo) localiza- 
tion radius defined by Eq. ( 13"7j) remains finite (but exponentially large) for any infinitely 
weak disorder, signalling the absence of Anderson transition in such system, similarly to 
the case of the usual single particle theory. However, as we shall discover below, local- 
ization radius in finite size systems diverges at some critical disorder, which is defined by 
the system size L. Qualitatively, this critical disorder is determined by the condition that 
localization radius of infinitely large system becomes comparable to characteristic sample 
size Rf^°° ~ L. Thus for finite two-dimensional systems Anderson transition in fact 
exists, as well as metallic phase for disorder below some critical value. In the following, 
under the term "correlated metal" phase we shall imply precisely such phase for finite 2D 
systems. 

Below we discuss most relevant results of DMFT+E approach for 2D Anderson- 
Hubbard model at half-filling on a square lattice with rectangular bare density of states 
with the bandwidth W = 2D: 



which corresponds right to a 2D case. 

Density of states and optical conductivity. 

Calculations show that qualitative behavior of the density of states in 2D is com- 
pletely analogous to that discussed above in three-dimensional case. Some quantitative 
distinctions are due to the different model of the "bare" density of states (132"]) . lead- 
ing in particular to larger (than in three-dimensional case) critical Hubbard interaction 
U C 2 ~ 1.83W, corresponding to Mott metal-insulator transition in absence of disorder 

5 These results for different temperatures were calculated by N.A. Kuleeva. 




s\ < D 
e\ > D ' 



(42) 
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and U& ~ 1A2W, which limiting below the region coexistence of metallic and insulating 
phases. Similar to three-dimensional model for U > U C 2 (i.e. for Mott insulator without 
disorder) increase of disorder leads to restoration of quasiparticle peak in the density of 
states. However, in this case such behavior does not signal in general the transition to a 
correlated metal state, as for infinitely large systems we are dealing here with correlated 
Anderson insulator. 

Optical conductivity behavior in a wide frequency range is also qualitatively is quite 
similar to that in three-dimensional model. But for infinite 2D model zero frequency 
conductivity always disappears (in the zero temperature limit) and, in contrast to d — 
3 case [27], even at very weak disorder the peak in optical conductivity lies at finite 
frequency. In the ladder approximation, which does not contain localization corrections, 
the usual Drude peak is observed at zero frequency and conductivity at u = is finite. 

In more detail results for DOS and optical conductivity in 2D model can be found in 
Ref. [69]. 

Localization radius and phase diagram of 2D Anderson-Hubbard model at 
half-filling. 

On the left side of Fig. [13]we show dependence of conductivity on disorder A for a finite 
but quite low frequency uj = 0.00005-D. Circles show results of ladder approximation, 
triangles - self-consistent theory of localization. Curve 3 qualitatively coinciding with 
ladder approximation was obtained from classical Drude expression: 

<r(u) = a (0)-£-, (43) 

where static conductivity cr(0) = e 2 N(0)D « y^-, N(0) - density of states at the 
Fermi level, D is Drude diffusion coefficient. Impurity scattering rate was taken as 
7 = 7rA(0)A 2 « 2I) A 2 . Significant contribution from localization corrections to conduc- 
tivity at finite frequency (noticeable distinction of curve 2 from 1 and 3) appears only 
when conductivity reaches values of the order of minimal metallic conductivity do = %■ 
(which is taken as unit of conductivity on figures). One should note that exactly in this 
range of disorder, as we shall see below, Anderson metal-insulator transition takes place 
(localization radius diverges) for 2D systems of reasonable finite sizes. 

Also in Fig. [13] we show dependences of localization radius logarithm following from 
( I37p (right scale) as function of disorder: curve 1 - infinite sample, curves 2 and 3 are 
for the finite size samples with L = 10 8 a and L = 10 5 a correspondingly. It is seen that 
localization radius grows exponentially as disorder diminishes and remains finite in the 
infinite 2D system, where Anderson transition is absent. For the finite systems localization 
radius diverges at critical disorder determined by system size, demonstrating the existence 
of an effective Anderson transition. As can be seen from Fig. [131 qualitatively the 
critical disorder is defined by condition that localization radius of infinite system becomes 
comparable with characteristic sample size R^ 00 ~ L. It should be noted that within our 
approach localization radius practically does not depend on U (in contrast to e.g. |61j), 
which leads to independence of critical disorder of correlation strength U in 2D finite 
size systems. Analogous situation is realized also for three-dimensional systems [27]. In 
general, it is of course a drawback of our approximations. 

On the same left part of Fig. [13] we plot the dependence of static conductivity on 
disorder strength in finite samples with sizes L = 10 8 a and L = 10 5 a (curves 4 and 
5 correspondingly). In the finite size systems with weak disorder static conductivity is 



27 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 




0.7 

0.6 

0.5 
0.4 

Q 

S °' 3 
0.2 

0.1 

0.0 



6=2 


U (A) ^ 


- Correlated Anderson Insulator 


/ 

if' VJ A U 


■ 










Mott Insulator 


■ "Correlated Metal" I <y 




j. 1 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
U/2D 



Figure 13: ON the left side - low finite frequency (to = 0.00005D) conductivity dependence 
on disorder strength A at U/2D = 1. Circles (curve 1) show results of ladder approximation, 
triangles (curve 2) - self-consistent theory of localization. Curve 3 (practically coinciding with 
ladder approximation) is obtained from Drude formula (|43H . Curves 4 and 5 - static conduc- 
tivity of finite samples with sizes L = 10 8 o and L = 10 5 a correspondingly. On the right scale 
localization radius logarithm versus disorder strength A: infinite size sample - curve 1, finite 
samples with sizes L = 10 8 a and L = 10 5 a - curves 2 and 3 |69| . 

On the right - phase diagram of 2D paramagnetic Anderson-Hubbard model at zero temperature 
[69] . Mott insulator region boundary f/ C 2(A) and coexistence region boundary U c i(A) were 
obtained from density of states behavior. Hatched part - region of an effective Anderson metal- 
insulator transition in finite size systems. 

not zero (metal) and gradually goes down with increase of disorder. It becomes zero at 
critical disorder, where localization radius diverges in the sample of corresponding size. 
Static conductivity of finite size samples within our approximation does not practically 
depend on correlation strength U. Significant difference of static conductivity from low 
finite frequency conductivity observed in Fig. [13] is related to exponential smallness of 
the frequency range with localization behavior of conductivity, mentioned above. 

Let us consider now the phase diagram of 2D paramagnetic Anderson- Hubbard model 
at half-filling, obtained from DMFT+S calculated densities of states and from the analysis 
of localization radius behavior in finite size 2D systems. Such phase diagram in coordinates 
disorder of A and correlation strength U is shown on the right part of Fig. [131 

Hatched stripe corresponds to the region of effective "metal"- Anderson insulator 
transition. Boundaries of this region are denned by divergence of localization radius in 
finite samples with characteristic sizes L = 10 5 a (upper boundary) and L = 10 8 a (lower 
boundary) (see left side of Fig. [13]). One should emphasize that further increase of the 
system size, e.g. ten times, up to L = 10 9 a, leads only to quite insignificant drop of 
critical disorder. In other words, it slightly shifts down low boundary of hatched stripe 
(on the right part of Fig. [T3]) - characteristic region of an effective Anderson transition 
in finite size systems. 

Curve ?7 c2 (A) computed from density of states behavior defines the boundary of Mott 
transition. Transition criteria is the disappearance of central quasiparticle peak in the 
density of states N(e) together with gap opening on the Fermi level. Similarly to three- 
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dimensional model, decrease of U starting from insulating phase leads to Mott transition 
U = U c i(A) < C/ c2 (A) and coexistence (hysteresis) region is observed on the phase dia- 
gram between curves U c i(A) and [/^(A) (Fig. [T3|) . In analogy with three-dimensional 
case we can guess that the ratio of Hubbard interaction and effective bandwidth ^ c /' c2 ff i ? 
controls Mott metal-insulator transition and is a universal constant, which does not de- 
pend on disorder and obtain qualitative dependencies U* 2 (A) for 2D model, which is 
plotted by dotted curve on Fig. [I3J It is seen that in contrast to d = 3 case [27] the 
dependence of U C 2(A) obtained from straightforward calculations of densities of states 
significantly differs from the qualitative U* 2 (A) dependence. Apparently it is related to 
important change of density of states lineshape (at U — 0) as disorder A grows, which is 
absent for semielliptic band in d = 3 case. 

3.3 Singularities of electron dispersion in strongly correlated 
systems in DMFT and DMFT+E approaches. 

3.3.1 Cusps ("kinks") in electron spectra. 

Electronic properties of crystalline solids are determined by single-particle and collective 
excitations of electron subsystem and their interactions with each other. These excitations 
are characterized by energy E and momentum k which are related to each other via 
dispersion (spectrum) E^. Interaction between single-particle and collective modes can 
result in noticeable bends (cusps) of dispersion E^ - the so called "kinks" . The lineshape 
and energy position of these features characterize interactions in the manybody system. 
For example, kinks in electron dispersions experimentally observed by ARPES in copper 
oxides, with energies 40-70 meV below the Fermi level, are evidence of electron-phonon 
[ZSl ESI [SHI ED E21 E31 El] or electron-magnon (spin-fluctuation) interactions [851185]. 

At the same time ARPES experiments detect kinks in electron dispersion for a number 
of different systems at essentially higher energies (up to 800 meV) jSTJ EH EH] . Physical 
nature of these kinks remains unknown. Among other attempts to find the explanation of 
these electron spectra anomalies a novel purely electronic mechanism of kinks formation 
was proposed in Ref. [90]. This mechanism is applicable to strongly correlated metals, 
where spectral function contains well developed Hubbard subbands, together with central 
quasiparticle peak as, for example, in transition metal oxides. Energy location of these 
kinks apparently determines the range of applicability of Landau Fermi liquid theory. 

As is well known [HI] in general case interaction results in finite life time of excitations 
in the system, so that E^ becomes a complex function. For electron systems with Coulomb 
interaction Landau Fermi liquid theory proves the existence of weakly damped fermionic 
quasiparticles for low enough temperature and in the narrow energy interval around the 
Fermi surface [91]. Beyond the Fermi liquid regime the concept of quasiparticles with 
well defined dispersion, strictly speaking, is not applicable since quasiparticle lifetime is 
too small. However, in recent years ARPES experiments convincingly show the existence 
of essentially k-dependent (though with rather broad maxima) single-particle spectral 
function behavior, pretty far from the Fermi level, despite the fact that st these energies 
one can not speak about well defined quasiparticles. In this case we understand as particle 
dispersion precisely this k-dependence of spectral function maxima which replaces the 
usual notion of quasiparticle spectrum. 
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3.3.2 Kinks of purely electronic nature 

Let us consider in more details the new mechanism of kinks formation in electronic dis- 
persion of strongly correlated systems, which does not assume an interaction of electrons 
with phonons or other excitations [90]. In order to understand the nature of this mecha- 
nism at the beginning we shall examine weakly correlated system described by standard 
Fermi liquid theory. Because of large enough quasiparticles lifetime close to the Fermi 
level and weakness of correlations, in the first approximation interaction leads to a sim- 
ple renormalization of initial dispersion of (noninteracting) quasiparticles e^, which is 
replaced by = Zp^e^, where Z-fl is electron Fermi liquid mass renormalization coeffi- 
cient. It changes the slope of dispersion in the vicinity of the Fermi level, compared to a 
bare one. However, if we consider an electron with energy far away from the Fermi level, 
then in case of weak interaction one can expect that its dispersion practically does not 
change: w e^, though the damping here can be pretty large. In this sense one can 
say that electron-electron interaction by it self can result in formation of bends (kinks) 
in the generalized dispersion (of the spectral function), and position of such kink on the 
energy scale is defined right as an energy, where simple Fermi liquid picture becomes 
inapplicable. However, in weakly correlated metals Zfl < 1, so that the slope of E^ will 
be changed insignificantly, making such kinks hard to observe. 

This picture can be essentially different in strongly correlated systems, where the 
value of Zfl can be much less than one, thus making kinks more pronounced. Strong 
correlation interactions give rise to strong spectral weight redistribution within the single- 
particle spectral function due to formation of Hubbard subbands. Moreover, as we know 
in strongly correlated metals there is also quasiparticle peak which appears close to the 
Fermi level, between Hubbard subbands. Below it will be shown that "usual" Fermi liq- 
uid quasiparticles exist in such systems only in a very narrow vicinity of the Fermi level, 
while beyond the Fermi liquid regime (but still "inside" of quasiparticle peak) exists some 
intermediate regime with dispersion E^ pa Zcpek, where renormalization factor Zqp is 
defined by central peak spectral weight and its value essentially differs from Zfl- At 
these intermediate energies, which are much smaller than interaction energy, electrons or 
holes are strongly interacting and their dispersion differs from both Fermi liquid one and 
noninteracting one. In this energy range one can speak about "intermediately" correlated 
situation, when Z FL < Z CP < 1. Consequently, at some energies ±u+ inside the quasi- 
particle peak transition from renormalization Zfl to Zcp takes place. That leads to an 
emergence of kink in electronic dispersion and the energy position of this kink is directly 
related to limits of applicability of the usual Fermi liquid theory. Lets emphasize that this 
mechanism gives rise to kink formation without any additional interaction with phonons 
or other collective modes. The only necessary condition for such kinks formation is the 
presence of strong electron-electron correlations in the system. 

For microscopic description of electronic kinks lets consider Hubbard model, which 
will be analyzed in the framework of standard DMFT, using numerical renormalization 
group (NRG) to solve effective Anderson impurity problem. For simplicity we consider 
single-band Hubbard model at half-filling. Strongly correlated regime in Hubbard model 
occurs when interaction value becomes of the order of bare bandwidth U ~ W. Con- 
sider as an example computational results shown in Fig. [T3J It is clearly visible that 
dispersion undergoes from Fermi liquid regime (line 1 on Fig. HI]) into described above 
"intermediate" regime (line 2 in Fig. H%|) with formation of well defined cusps in dispersion 
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Figure 14: Cusps (kinks) in electronic dispersion E'k for the case of strongly 
correlated systems. On the contour plot intensity map of DMFT spectral function 
A(k, oj) for the Hubbard model on a cubic lattice is presented. Hubbard interaction value 
is £7=3.5 eV, bandwidth W ~ 3.64 eV, band filling n—1, calculated value of Fermi liquid 
renormalization is Zfl=0.086, temperature is T = 5 K. Close to Fermi level maxima of 
A(k,oj) (white dots) correspond to renormalized dispersion Ey_ = Zp^e^ (line 1). At ener- 
gies \u\ > u+ spectral function A(k,u) keeps its shape but with different renormalization 
E k = Z CP e k - c sgn(£ k ) (line 2). Values of w*=0.03 eV, Z CP = 0.135, and c = 0.018 eV 
are calculated in [90] from the values of Zfl and corresponding to black line. Inset on 
the right shows in details part of dispersion in the T-R direction marked out with white 
rectangle, kinks (cusps in dispersion) at energies ±u+ are pointed by arrows. Black lines 
- bare (noninteracting) quasiparticle spectra. 



(kinks) at energies = 0.03 eV. In some of the high symmetry directions (e.g. around 
X-point, Fig. [T4"l) dispersion has quite small slope close to the Fermi level making kinks 
less pronounced. 

Let us have a look on this situation from the point of view of single-particle Green's 
function behavior. From general considerations it is clear that any cusps of dispersion 
caused by interaction appear because of corresponding behavior of the self-energy or, 
more precisely, that of Re£k(a;). In the majority of real physical systems k-dependence 
of self-energy, apparently, is less important in comparison with w-dependence and thus the 
neglect of self-energy k-dependence Sk(w) = S(w) is, more or less, good approximation. 
In the framework of DMFT this statement is exact. Self-consistent expression for self- 
energy of Hubbard model within DMFT can be written as: 

E(w)=w + ^-l/G(w)-A(G(w)), (44) 

where G(cu) = 4 G(k, uj) is the local (averaged over k) Green's function, A(G) is 
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frequency dependent hybridization function expressed via G{uS) %. Hybridization function 
describes quantum-mechanical coupling between an electron on a given site and other 
sites of the system. 

In Fig.[15]we show frequency dependence of spectral function (DOS) A(u) = — lmG(u))/ir 
calculated for the same model parameters as in Fig. and demonstrating a typical three- 
peak structure. Corresponding real parts of Green's function G(u>) and self-energy S(o;) 
are shown in Fig. [T5b and Fig. [T5b . 

Kinks in Re£(co>) are directly related to three-peak structure of integrated spectral 
function (DOS) A(u) (its minima are at energies ±f2) and define some new quite small 
energy scale. To this end Re[G(u;)] should have maxima and minima in the energy interval 
±w max , i.e. "inside" the central peak (Fig. [15b) . In its turn it directly leads to kinks 
formation in ReE(u;). Self-energy £(a>) consists of two contributions: oj + \i — l/G(u) 
and —A(G(u)). Expression Re[u; + fi — 1/G(u)} is linear within quite large energy range 
\oj\ < Q (Fig. [T5H) . while the value of —Ke[A(G(uj))) is proportional to —Re[G(u)], at 
least in the first order of corresponding momenta expansion, only in the small energy 
interval \u\ < u max . The sum of these two contributions gives rise to cusps in real part 
of self-energy at energies ±u;*, where u+ = (y/2 — l)u; max . At this energy Re[G(o;)] has 
maximal curvature (these points are marked out by circles in Fig. [T5b). Thus the Fermi 
liquid regime, when the slope of the real part of self-energy is described by <9Re£(w) / du = 
1 — 1/Zfl, exists only within a narrow part of the central peak, namely in the energy 
interval \u\ < w*. At higher (intermediate) energies the slope of ReS(w) will be given 
by dKeTj(u)/du = 1 — 1/Zqp- As a result, effective dispersion will manifest kinks at 
energies u*. 

This analysis also helps to understand why outside Fermi liquid region demon- 
strates another value of renormalization of noninteracting electron dispersion, given by 
Zqp with small offset c. This behavior is determined by the value of main contribution 
to self-energy oj + fi — 1/G(u) at energies u+ < \oj\ < Q i.e. within the central peak of 
DOS. Values of u*, Z CP , and c can be expressed via Z FL together with characteristics of 
noninteracting electron density of states. One can show that = Zfl (y/2 - 1)D, where 
D is the halfwidth of the bare band (details see in work [90]). If correlations are weak and 
Zfl < 1 kinks positions in E^ practically coincide with the edges of bare electron band, 
which makes them almost unobservable. On the other hand, in the strongly correlated 
regime (Zfl ^ 1) kinks energy u*/D oc Zp^ approaches the Fermi level inside the central 
peak, which width diminishes with the increase of correlations as Q/D oc V^fl [94] . 

For the first time, these purely electronic kinks were observed in LDA+DMFT calcu- 
lations for SrV03 system [92] • Definition of energy scale u± contains only parameters of 
initial band structure, which can be obtained (for realistic systems) via band structure cal- 
culations, together with Fermi liquid mass renormalization Zfl = 1/(1 — <9ReX(0)/<%j) = 
m/m*, which can be experimentally determined from specific heat or spin susceptibility 
measurements In particular, in Ref. [93] it was shown that kinks of electronic nature can 
cause corresponding cusps in the linear (in temperature) term of specific heat of strongly 
correlated metals, which was, apparently, observed in LIV2O4 heavy fermion system. 

6 In DMFT A(G) is defined by G(w) = G q {lo + fj, - E(w)), i.e. G (A(G) + 1/G) = G, where G (lu) is 
noninteracting local Green's function. 
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Figure 15: Green's function characteristic for strongly correlated system with parameters 
given in the caption of Fig. HU a - three-peak structure of spectral function A(u) = 
—lmG(u)/iT with minima at energies ±f2 = 0.45 eV. b - real part of Green's function 
with relative minima and maxima at energies ±a; max located "inside" quasiparticle peak, 
c - real part of self-energy with cusps at energies ±u± (circles) situated in the points of 
maximal curvature of ReG(u) (u± = 0.4w max = 0.03 eV). d - real part of uj — 1/G(u) 
(solid line) contributing to self-energy which is linear, in general, in the energy interval 
\co\ < Q. Change of energy dependence of — Re[G(u;)] at energies ±u+ defines the location 
of kinks (cusps) in electronic dispersion. 
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3.3.3 Role of electron-phonon interaction 

The previous section material inevitably brings us to the question about the relation and 
mutual influence of kinks of electronic nature and the "usual" kinks in the electron disper- 
sion induced by electron-phonon interaction. This is related to a more general problem 
of joint influence between strong electronic correlations and electron-phonon interaction. 
In fact, the history of such investigations is relatively long and one of the most popular 
models of electron-phonon interaction in strongly correlated systems is Hubbard-Holstein 
model. Hubbard model describes local Coulomb interaction on a lattice. On the other 
hand, Holstein model describes linear interaction of conduction electrons with local (Ein- 
stein) phonon modes [95]. Studies of Hubbard-Holstein model were performed in the 
framework of conventional DMFT [5], in particular with the use of numerical renormal- 
ization group (NRG) [15] as "impurity solver". Reducing of Hubbard-Holstein model 
to Anderson-Holstein impurity problem was first performed by Hewson and Mayer [96] . 
They showed that using NRG one can calculate the total electron-phonon contribution to 
self-energy, thus achieving nonperturbative solution of Hubbard-Holstein model not only 
with respect to Hubbard interaction, but also with respect to electron-phonon interaction. 
Let us note that the general structure of DMFT equations in this approach is preserved. 

However, until recently there were no studies of strongly correlated electrons interact- 
ing with Debye phonons. It is quite surprising in view of rather wide discussion of the 
physics of kinks in electron dispersion observed in ARPES experiments in high-T c super- 
conducting oxides (El)] ■ The origin of these kinks is typically attributed to electron-phonon 
interaction [ST]. Problem of kinks formation in electron dispersion caused by electron- 
phonon interaction in strongly correlated systems was shortly discussed in the framework 
of Hubbard-Holstein model Ref. [971 EH]- bn this section we overview DMFT+S results 
for Hubbard model with added interaction with Debye phonons, under the assumption of 
the validity of Migdal theorem (adiabatic approximation). This approximation is reason- 
able for electron-phonon coupling constant A < Epjup, ~ 10, where Ep is Fermi energy, 
ujd is Debye frequency. 

To consider electron-phonon interaction in the strongly correlated system we introduce 
self-energy £k( £ ) — k), appearing in the usual Froehlich model. According to 

Migdal theorem, in adiabatic approximation we can restrict ourselves to the simplest 
first order contribution to S p / l (e,k). Main advantage of adiabatic approximation is the 
possibility to neglect any vertex corrections from electron-phonon interaction, which are 
small over adiabatic parameter ^ <C 1 [99] . 

In fact, Ep/^e, k) in this approximation has only weak momentum dependence, which 
can be neglected and we have to account only for essential frequency dependence. Direct 
calculations (see e.g. [29]) in the case of Debye phonon spectra wo(k) = w|k| for |k| < 
where u is sound velocity, give: 

—■in 2 r + °° rli i , ,2 . ,2 

Z ph (e) = -^ ^{ujl + unn\^^\+i^6{ul-^)}l{e + u) 1 (45) 

4d>c J-oo 2tt u 2 

where g is the usual electron-phonon interaction constant /(e) = d ^ ^ , E £ = 
e— S(e) — T.phi^e) and u c = ppu is characteristic frequency of the order of Debye one. In case 
of semielliptic bare DOS Nq(e) with halfwidth D one obtains 1(e) = -j=p(E £ — \J E 2 — D 2 ). 
It is convenient to introduce dimensionless electron-phonon coupling constant, which for 
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Figure 16: Comparison of DMFT (dashed line) and DMFT+E P /, (solid line) densities of states. 
On the left side - for strong (upper panel, U /2D =1.25) and weak (lower panel, U/2D =0.625) 
Hubbard interactions. On the right side - densities of states evolution close to metal- insulator 
transition. Dimensionless electron-phonon coupling constant is A=0.8 |100|, 1101] . 



this model can be written as [29] : 

X = 9 2 No(e F )^- 2 . (46) 

To simplify calculations we neglect renormalization of phonons because of electron-phonon 
interaction [29], assuming in the following that phonon spectrum is fixed by experiment. 



3.3.4 Electronic and phonon kinks within the DMFT+S approach 

Let us focus on most interesting DMFT+S results in this model, referring the reader for 
details to Refs. [1001 1101] . Here we present results obtained for the case of interaction of 
electrons with Debye phonons (results for Einstein phonons is analogous [101j ). 

Comparison of DMFT and DMFT+Eph densities of states with electron-phonon inter- 
action for strong (U/2D=1.25) and weak (U/2D=0.625) Hubbard interaction is presented 
on upper and lower panels on the left side of Fig. [16j Dimensionless constant P6|) ex- 
ploited in these calculations was chosen to be A=0.8 and Debye frequency ajD=0.125D. 
In both cases we observe some spectral weigh transfer caused by electron-phonon interac- 
tion. For U/2D=1.25 (upper panel of Fig. [TH]) we see well developed three-peak structure 
typical for strongly correlated systems. In the energy interval ±u/j around the Fermi 
level (which is zero energy in all figures) practically there is no difference in the shape 
of quasiparticle peak in the DOS obtained within DMFT and DMFT+S p /j. But beyond 
this interval DMFT+Ep^ quasiparticle peak becomes significantly wider because of partial 
transfer of spectral weight from Hubbard bands. This broadening of quasiparticle peak 
in DMFT+Ep/! leads to a delay of metal-insulator transition as we shall see later. 

In case of U/2D=0.625 no clearly distinguishable Hubbard bands are formed and we 
observe only some side wings in the DOS. Redistribution of spectral weight on the lower 
panel of Fig. UHlis not very dramatic despite qualitative difference with the U/2D=1.25 
case. Main distinction between DMFT and DMFT+E P /, results occurs here in the interval 
±ujd, where formation of a "cap" in DOS is observed, due to electron-phonon interaction. 
Corresponding spectral weight goes to the energies around ±U where Hubbard bands 
start to form. 
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On the right panel of Fig.. US] we compare behavior of densities of states in DMFT 
without phonons and in DMFT+Ep^, for different values of U/2D parameter nearby 
Mott-Hubbard transition. At U/2D=1.56 both standard DMFT and DMFT+£ p/l give 
insulating solution. However, there are some distinctions between these solutions. In 
DMFT+Eph Hubbard bands are lower and wider than in DMFT since additional (electron- 
phonon) interaction is included. With decrease of U at U/2D=1.51 and 1.47 we observe 
that DMFT+Sp/j results correspond to metallic state (with narrow quasiparticle peak on 
the Fermi level), while DMFT without phonons still shows insulating solution. Only about 
U/2D=1.43 both DMFT and DMFT+S p/l results for DOS correspond to the metallic 
state. Thus, under increase of U the finite value of electron-phonon interaction slightly 
delays Mott-Hubbard transition from metallic to insulating phase. This result is analogous 
to one obtained within the Hubbard-Holstein model for weak electron-phonon coupling 
[T02lfT03lfT03] . 

Therefore, moderate electron-phonon interaction brings to rather insignificant changes 
of electron density of states both in correlated metal and in Mott insulating state, only 
slightly delaying transition from metal to insulator with the growth of U. 

Now we turn to the source of sharp slope changes of electron dispersion (kinks). It is 
well known that in general case kinks are formed because of interaction of electrons with 
bosonic modes. In case of electron-phonon interaction, typical energy of the kink is about 
Debye (or Einstein) frequency. Above we have shown that in strongly correlated metal 
kinks of purely electronic nature can arise [90]. Energy of such kink for semielliptical 
bare DOS is u* = Zfl(V^ — 1)-D, where D is the halfwidth of the "bare" band and 
Zfl — (1 — 9B Qe \e=E ) _1 * s Fermi liquid renormalization factor. Roughly speaking co* is 
defined by halfwidth of the quasiparticle peak in DOS. 

Kink of electronic nature is quite smooth and its observation is rather difficult. DMFT+E p /j 
calculations show that electronic kinks are hardly detectable on the background of phonon 
kinks and fine "tuning" of model parameters is necessary pick them out. First of all, it 
is necessary to guarantee that cup <C u* (in other cases smooth electronic kinks will be 
practically indiscernible against kinks from electron-phonon interaction). For U/2D=1 
and U=3.5 eV we have u* ~ 0.1D, while Debye frequency can be taken quite small e.g. 
u>d ~ 0.01D. In order to make phonon kink pronounced enough at such relatively low 
Debye frequency one needs to increase electron-phonon coupling constant up to A=2.0. 

To demonstrate the possibility of coexistence of both types of kinks in the spectra, let 
us consider the energy dispersion for simple cubic lattice with nearest neighbors hopping 
only. Most convenient is to discuss the high symmetry direction T — (tt, tt, tt) of the 
Brillouin zone [HD]- On the left panel of Fig. [TT] we show electronic dispersion along 
this direction close to the Fermi level. Line with diamonds is the electron spectrum 
of standard DMFT without phonons. Lines with circles presents DFMT+Ep^ results. 
Electronic and phonon kinks are marked by arrows. In general case, kinks from electron- 
phonon interaction dominate for most typical model parameters making electronic kink 
observation predicted in [HD] quite difficult. 

In conclusion we give the picture of phonon kinks evolution depending on the value 
of Hubbard interaction U. With decrease of U/2D ratio Fermi velocity goes down and 
kink position on the momentum axis shifts farther away from pp, while the kink energy 
remains about ujd- This behavior follows from direct DMFT+E p /j calculations |1QQ[ 1101] 
and is shown on the right panel of Fig. [T7J In the case of interaction with Einstein phonons 
results are quite analogous [101J. 
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Figure 17: Quasiparticle dispersion with kinks around the Fermi level. On the left side - 
dispersion along T — (ir, it, it) high symmetry direction in the Brillouin zone for the case of 
simple cubic lattice and bare dispersion with nearest neighbors hopping only: DMFT (curve 
with diamonds) and DMFT+Ep/j (curve with circles) (A = 2, ljd = 0.01D). On the right side 
- DMFT+Sp/j phonon kinks evolution for different values of Hubbard interaction U/2D = 0.5, 
0.75, 1.0; A = 0.8, oj d = 0.1D [TOO] . 

4 Electronic structure of real strongly correlated sys- 
tems: LDA+DMFT and LDA+DMFT+E 

At present the most advanced ah initio (i.e. ideally without any fitting parameters) com- 
putational scheme of electron spectra of realistic strongly correlated systems is LDA+DMFT 
method [21 [16]. LDA band structure in this method is used to obtain "noninteracting" 
starting Hamiltonian, while strong correlations are accounted afterwards within DMFT. 
In fact, LDA+DMFT computational scheme combines two scientific areas: "realistic" 
band structure calculations and traditional model approaches, which were essentially sep- 
arated from each other before. Without DMFT just in LDA there is no way to describe 
strongly correlated systems, while without LDA calculations manybody methods can not 
be material specific. Below we briefly discuss the basics of LDA+DMFT and generalized 
LDA+DMFT+S methods. 



4.1 Density functional theory (DFT). Local density approxima- 
tion (LDA). 

In Born-Oppenheimer adiabatic approximation |105j . neglecting relativistic effects, elec- 
tronic properties in solid state physics are described by the Hamiltonian: 



H = J2fd 3 r ^ + (r, 



cr 



ll 2 



2m, 



-A + V^ on (r) 



+ \^ J " d3rd3r ' ^ + (r,<x)^+(iV) Ke(r-r') *(r', a')$(r,a). (47) 
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Here ty + (r,o~) and 4?(r,a) are creation and annihilation operators of electron with coor- 
dinate r and spin a, A is the Laplace operator, m e is electron mass, e is electron charge, 

vUr) = -e 2 J2r\i > 7 -M = 7Eri (48) 

/ — ' r — Rj 2 r — r 

i r^r' 1 

denote single particle potential created by all ions i with charge eZi located at given 
positions Rj and electron-electron interaction. 

Although the "first principle" Hamiltonian (T4T]) is easy to write down, it is impossible 
to solve corresponding quantum mechanical problem exactly. This is the reason to make 
substantial physical approximations. In particular, density functional theory (DFT) is 
based on Hohenberg-Kohn theorem [106] (see also the review |107j ). which asserts that 
ground state energy is the unique functional of electron charge density, which is minimal 
for the equilibrium electron density in the ground state: 

E [p] = E kin [p] + E ion [p] + £ Hartr c C [p] + E xc [p] , (49) 

where the Hartree energy -EHartrec[p] = §/ d 3 r' d 3 r V^ e (r — r') p(r')p(r) and potential 
energy of ions Ei on [p] = J d 3 r Vi on (r) p(r), can be directly expressed via electron charge 
density. The term -Ekirjp] denotes the kinetic energy of electrons and E xc [p] is unknown, in 
general, exchange-correlation term, containing electron-electron interaction energy beyond 
the Hartree term. In fact, all peculiarities of the manybody problem are transferred into 
the computation of E xc [p\. 

In practice, instead of minimization of E[p] over p, minimization is performed usually 
over some set of orthonormal functions tpi, related to p via expression: 

N 

p(r) = 5>,(r)| 2 . (50) 
i=i 

Introducing arbitrary Lagrange parameters £j and requiring 

5 



j£[p]+ £i [l- J d 3 r\^(r) 



one gets Kohn-Shem equations [108J: 



-A + Mon(r) + J d 3 r'V cc {v-v')p{v') 



(51) 



V? t (r) = £i <pi(r). (52) 



2m e J ' <^p( r ) 

Formally this equation coincides with single particle stationary Schroedinger equation. 
Electron kinetic energy corresponding to charge density of the ground state is given now 
by the expression 

JV 

£kin[Pmin] = ~ ^(^1 ^ 2 A/ (2m e )\(pj) , (53) 

1=1 

where ifi are self-consistent (spin degenerate) solutions of equations (152]) and ( |50l) . corre- 
sponding to the smallest single particle energy |109j . 

Most common way to calculate E xc [p] is to use the local density approximation (LDA). 
It approximates the functional E xc [p] by the local charge density functional: 

E xc [p]^ J rf 3 r/)(r)i DA Wr)) (54) 
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Explicit expression for e^ A (p(r)) can be found in the framework of different models e.g. 
from numerical analysis of "jellium" model (electronic gas on a positive ionic background) 

mm- 

In fact, LDA approximation corresponds to the replacement of the Hamiltonian ( 1471) 

by 



#LDA = Yl J^ T 



(J 



h 2 



+ 



2m e 
5E^[p} 



A + V ion (r)+ / dVp(r')^ee(r-r' 



8p(r) 



(55) 



In practical calculations electron field operators are expanded over some set of atomic- 
like functions (basis set) &u m {i denotes lattice site, I is orbital quantum number, m is 
magnetic quantum number). In this representation 



*+(r,a) = £cf/ m <Wr) 



(56) 



ilm 



and the Hamiltonian fl55l) is rewritten as 



H 



LDA 



^ ] (fiilm,jl'm' £ ilm ^ilm + t'ilmjl'm' ^ilm^jl'm') • 
ilm, j I'm' ,u 



(57) 



Here ff ilm = c[l m c a ilm is electron density operator on a given orbital, while matrix elements: 



t 



ilm,jl'm' 



ilm 



h 2 A 

2m P 



+ V Um (r) + j dVp(r')Ke(r-r') 



8p(r) 



Qjl'm' 



(5f 



in case of ilm ^ jl'mf define effective hopping integrals and single particle energies e# m 
are given by corresponding diagonal expressions in case of identical indices. On this step 
purely analytic work ends and numerical calculations follow within the chosen basis set, 
e.g. linearized muffin-tin orbitals (LMTO). Specific expressions for matrix elements f )58|) 
within the LMTO basis are presented in Ref. [11 lj . 



4.2 LDA+DMFT computational scheme 

The primary importance for strongly correlated materials is the onsite Coulomb interac- 
tion between d- or /-electrons, since this contribution to interaction energy is the largest 
one. The largest nonlocal contribution is density- density type interaction between nearest 
neighbors, where the main contribution comes from Hartree term (see [112] and |113j ). 
which is already taken into account in the LDA. Moreover in LDA some part of exchange- 
correlation interactions is taken into account via effective potential SE ^ r ^ ■ 

In order to take into account strong local Coulomb interaction one should supplement 
Hamiltonian f )57|) with approximate Coulomb matrix with most important parameters 
only [16J: U - onsite intraband Coulomb repulsion, J - exchange interaction and in- 
terorbital Coulomb repulsion U', acting on different electronic orbitals of the same site id, 



39 



where sits an atom with partially filled <i-shell (I - orbital quantum number, m - magnetic 
quantum number, a - spin index): 



m i=i d) l=l d 

(59) 

i=id,l=ld m^m' era' 
i=id,l=ld ma 

where for simplicity only density-density type interactions are left, and so called Kanamori 
parametrization is applied, when for the same orbitals (m = w!) direct Coulomb interac- 
tion is taken as U, while for different orbitals (in ^ m') this interaction is equal to U' . 
Because of rotational invariance of the single atom problem U' — U — 2 J and exchange 
interaction parameter does not depend on orbital index and is equal to J. 

Moreover, in (160]) the last Ae^ - term is added (the so called double counting correc- 
tion), which should correct for double counting, as some part of local Coulomb interaction 
is already included into Hlda- General expression for Ae^ via U and p is unknown. How- 
ever, there are several qualitative recipes to determine the value of Ae^, and which are 
employed in different modern LDA+DMFT calculations (detailed discussion can be found 
in Refs. |114[ 1115] ). The simplest physical assumption is that within the DFT Coulomb 
interaction energy can be written as: 

E DFT = -Un d {n d -1), (60) 

where n d is the total number of electrons on <i-shell and U is the average Coulomb inter- 
action (here we assume averaging over all orbital pairs ma, m'a' on a given site). Thus 
Aed is taken as: 

The values of interaction parameters U, J and U' can be obtained from the averaged 
Coulomb interaction U and Hund exchange parameter J. The averaged interaction U is 
related to U and U' parameters via the following relation: 

,-, u + (N„ b - 1)U' + (N„ b - 1)(U' - J) 

U = 2iV^l ■ (62) 

where N or b is the number of interacting orbitals. Since U and U' are not independent 
parameters U and J are sufficient to determine U [20j 1116] . 

For microscopic calculations of averaged Coulomb interaction different methods were 
developed, such as "constrained LDA" [117] or "constrained RPA" [1181 1119] . General- 
izations to calculate Hund exchange parameter also exist. Unfortunately, there are rather 
large discrepancies between the values of parameters obtained with these methods. The 
value of U appears to be strongly dependent on the basis set used (e.g. in the problem 
of screening of long-range part of Coulomb interaction). It is clear that introduction of 
all these essentially model parameters takes us quite far away from "first principle" ideal, 
though it is the best one can do at the moment to calculate band structure of solids with 
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transition metal atoms, where electron-electron interactions play the crucial role. In that 
sense it is probably more correct to speak about "modelling" of electron structure of such 
systems. 

Matrix elements of "noninteracting" Hamiltonian in the reciprocal space H® DA (k) can 
be calculated numerically at every point of the Brillouin zone, then the integrals over the 
Brillouin zone usually are calculated with tetrahedron method |120] . In case of relatively 
simple band dispersions, when analytical expression for H^ DA (k) dependence on k can be 
written explicitly, the values of hopping integrals can be found from LDA by projecting 
on corresponding Wannier functions [1211 I122j . Matrix elements of this Hamiltonian, 
i.e. single particle LDA energies without local Coulomb interaction can be written in a 
following way: 

(.H® DA (k)) q i mjq i[i m r = {Hjj D xi} S ))qlm,q'l'm' ~ ^qlm,q'l'm^ql,q d l d ^d n d- (63) 

where q is index of an atom in the primitive unit cell. 

The essence of the next step is to use in DMFT or DMFT+S loop (see section l2~7Lj) the 
local lattice Green's function (jSJ), determined by momentum integrated Dyson's equation 
of the form: 

Gqlm,q'l'm'{' jJ ) = TJ~ f dk [ U 5 q l m ,q'l'rn' ~ (-^LDA ()^))qlm,q'l'm' 
VB J 

+ bqlmh ^qlm,q'l'm'{u)]~ l , (64) 

where denotes the inverse matrix with indices n(=qlm), n'(=q'l'm'), while the 

integration is performed over the Brillouin zone of the volume Vb- 

Significant simplification of computations is achieved for the case of cubic lattice sym- 
metry, when the crystal field strongly splits <i-orbitals into the three-fold degenerate t2 9 - 
states and two-fold degenerate e 5 -states. In this special case both Green's function and 
self-energy become diagonal over orbital and spin indices. Then the calculation of local 
Green's function of the lattice problem can be performed as energy integration with the 
use of unperturbed density of states, which allows to avoid tedious integration over the 
Brillouin zone in (|64p and write: 

G(u) = G\u " = / de ^ ■ (65) 

In this case double counting correction Ae^ reduces to immaterial shift of the chemical 
potential and its particular mathematical form is irrelevant. 

4.3 Examples of LDA+DMFT calculations. 
4.3.1 Cubic perovskites CaV0 3 and SrV0 3 

In this section we consider examples of some LDA+DMFT calculations of electronic band 
structure of realistic compounds with strong enough electronic correlations. Transition 
metal oxides are ideal testing area to study electronic correlations in solids. Among these 
materials cubic perovskites have simplest crystal structure and thus can be viewed as 
a starting point to understand electronic properties of more complex systems. Usually 
3d-states in such materials form comparatively narrow bands of the width iy~2 — 3 eV, 
leading to strong electron-electron correlations. 
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Modern stage of experimental investigations of spectral and transport properties of 
strongly correlated 3d 1 transition metal oxides started from the work of Fujimori et 
al. [123] . The authors, apparently for the first time, discovered strongly pronounced 
lower Hubbard band in photoemission spectra, which could not be explained by standard 
methods of band structure calculations. In many of earlier works |124] 1125] 1126] 1127] . 
devoted to the properties of the series of compounds Sri^Ca^VOs with different val- 
ues of x rather controversial results were reported. While thermodynamic characteristics 
(Sommerfeld coefficient, electric resistivity and magnetic susceptibility) appeared to be 
more or less x independent, spectroscopic measurements data rather strongly changed 
as system transformed from x = (SrVOs) to x — 1 (CaVOs). These data indicated a 
transition from strongly correlated metal (SrVOs) to practically ideal insulator (CaVOs), 
with concentration range x — Y 1 in Sri_a;Ca x V03 being the boundary of Mott-Hubbard 
transition. Analysis of this problem was performed using the high penetration depth pho- 
toemission experiments by Maiti et al. [128], and similar experiments with high resolution 
photoemission by Sekiyama et al. [129 J. In particular, in the last work it was shown that 

(1) surface preparation technique is very important (preferable is cleavage method) and 

(2) energy of X-ray beam should be big enough to provide penetration depth of several 
elementary cells. At the same time high instrumental resolution should be guaranteed 
(about 100 meV in the work |129j ). Such improvement of photoemission methods lead to 
observation of almost identical spectra for Sr(Ca)VC>3 [128, 1129] . demonstrating agreement 
of spectroscopic and thermodynamic measurements. Results of these experiments agree 
also with earlier Is X-ray absorption spectra (XAS) obtained by Inoue et al. |130j . which 
differ only for energies slightly below the Fermi level in contrast to BIS data |127] . In the 
framework of single band Hubbard model with neglect of orbital structure of 3d-shell of V, 
Rozenberg et al. |131] modelled Sr 1 _ :r .Ca a; V03 spectra obtained by high penetration depth 
photoemission |128] using adjustable parameters. Later in Ref. |129] it was demonstrated 
that data Ref. [128] contained quite significant surface contribution. 

Below we present results of LDA+DMFT(QMC) calculations, performed without any 
adjustable parameters, both for spectral function and density of states of cubic SrVC>3 
and orthorhombic CaVOs perovskites. According to these both systems are strongly 
correlated metals, which are quite far away from metal-insulator transition boundary. 
Despite significantly smaller V-O-V bond angle in CaVOs, photoemission spectra of both 
systems are very similar and their quasiparticle peaks are almost identical. The results 
obtained agreed very well with modern high resolution bulk sensitive photoemission data, 
mentioned above. In the spectral function of SrVOs, obtained from LDA+DMFT(QMC) 
calculations, kinks of purely electronic nature at about 200 meV were observed, and later 
these kinks were observed experimentally. 

Results of LDA+DMFT calculations First of all from LDA calculated band struc- 
ture we extract single electron Hamiltonian H^ jDA with subtracted average Coulomb inter- 
action (to avoid double counting) [IB] . Supplementing i^DA with local Coulomb interac- 
tion between electrons we obtain the Hamiltonian ( 159]) for the material of interest. Since 
symmetry of CaV0 3 is close to cubic one, it is possible to simplify the calculations and 
use integration with LDA density of states N°(e), instead of integration over the Brillouin 
zone. In the Hamiltonian (155]) local intraorbital and interorbital repulsions and exchange 
interactions are taken into account explicitly as U, U' and J. The values of these interac- 
tions for SrV03 were calculated by constrained LDA method |117] with e 5 -states included 
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Figure 18: On the left side: LDA+DMFT(QMC) spectra for SrV0 3 (solid line) and 
CaV0 3 (dashed line) at T=300K (inset: temperature influence on CaV0 3 spectrum line- 
shape). On the right side: comparison of calculated (without adjustable parameters) 
LDA+DMFT(QMC) spectra for SrV0 3 (solid line) and CaV0 3 (dashed line) with high 
resolution bulk sensitive photoemission data (SrV0 3 - circles; CaV0 3 - rectangles) |129j 
(left picture) and Is- XAS-spectra: (SrV0 3 - diamonds, Cao.gSr .iV0 3 - triangles) |13Uj 
(right figure). Horizontal line - experimental background. 



into screening [132] . Obtained value of averaged Coulomb interaction is U = 3.55 eV 
(U = U' for t2 9 orbitals [IT] 1116] ) and J = 1.0 eV. Intraorbital Coulomb repulsion U is 
fixed by rotational invariance U = U' + 2 J = 5.55 eV. For CaV0 3 U was not calculated, 
since standard procedure of calculation of Coulomb interaction parameters between two 
t2g electrons screened by e g states is not applicable for distorted crystal structure, where 
e g and t<i g orbital are not separated by symmetry. On the other hand it is known that 
changes of local Coulomb interaction are usually much smaller than changes in density of 
states, which as shown above are weakly dependent on bond angle V-O-V. It means that 
U for CaV0 3 should be practically the same as for SrV0 3 . Correspondingly, the values 
U = 3.55 eV and J = 1.0 eV were used for both SrV0 3 and CaV0 3 . These values agree 
with other band structure calculations for vanadium compounds |132j and experimental 
data [133] . 

Further calculations with Hamiltonian (159]) were performed in the framework of DMFT 
with quantum Monte-Carlo method (QMC) [13J as an impurity solver. In QMC the 
Green's function was obtained for the imaginary time and then continued on the real 
time (frequency) by maximum entropy method [134] . In LDA+DMFT(QMC) calculated 
spectra for SrV0 3 and CaV0 3 shown in Fig. [18] (on the left side) we observe manifestations 
of correlation effects, such as formation of lower Hubbard bands near —1.5 eV and upper 
Hubbard bands at about 2.5 eV with well developed quasiparticle peaks on the Fermi level. 
Thus both SrV0 3 and CaV0 3 are strongly correlated metals. Difference of bare bandwidth 
(about 4%) only leads to a small additional spectral weight transfer from quasiparticle 
peak to Hubbard bands and slight changes of Hubbard bands positions. Obviously, both 
systems are close to Mott-Hubbard metal insulator transition. Manybody densities of 
states for both systems (Fig. [TS]) are similar but not identical. Indeed, SrV0 3 is a bit less 
correlated than CaV0 3 in agreement with difference of the LDA bandwidths. The inset 
in Fig. [18] shows that temperature influence on the spectrum is small for T < 700 K. 

In the middle part of Fig. \M LDA+DMFT(QMC) densities of states (obtained at 
T=300K, multiplied with the Fermi function at 20 K and broadened with Gaussian of the 



43 



width 0.1 eV to mimic experimental resolution [129J ) are compared with experimental 
photoemission data after subtraction of oxygen and surface contributions. In particular, 
the height and the width of these spectra are almost the same in SrVC>3 and CaVOs (with 
a bit of difference above the Fermi level). On the other hand, positions of lower Hubbard 
band differs quite noticeably. This distinctions might occur because of subtraction of 
(estimated) oxygen contribution, which can delete part of 3d spectral weight below —2 eV, 
as well as with uncertainties of U calculations. 

On the right side of Fig. [151 we show comparison with XAS data. Finite lifetime effects 
for holes are taken into account by broadening of theoretical spectra with Lorentzian of 
the width of 0.2 eV |135j . multiplication by inverse Fermi function (at T=80 K) and 
further broadening with Gaussian for experimental resolution 0.36 eV. General agreement 
of weights and positions of quasiparticle and upper Hubbard bands for t 2 g-band is good, 
including tendencies associated with transition from .SrV03 to CaVOs (in the experiment 
Cao.9Sro.1VO3). For CaV03 quasiparticle spectral peak weight is a bit smaller than in 
the experiment. In contrast to single band Hubbard model calculations LDA+DMFT 
accounts for peculiarities of the systems and reproduce strong asymmetry of the spectra 
close to the Fermi energy, including relative weights and bandwidths. These results give 
different interpretation of XAS as compared with Ref. |130j . where maximum at 2.5 
eV was associated with e g band and not with upper Hubbard band of t 2g band. Small 
differences of qusiparticle peaks (see Fig. H~8l) lead to different values of effective masses: 
m*/mo = 2.1 for SrVOs and m* /mQ = 2A for CaV03. These theoretical values agree with 
m* /mo = 2 — 3 for SrVOs and CaV03 obtained from de Haas- van Alp hen experiments 
and thermodynamic data [1241 H251 11261 H36j . Note that the effective mass for CaVOs 
determined from optical experiments is slightly larger: m*/mo = 3.9 |133j . 

4.3.2 Kinks in spectral function of SrVOs. 

Let us consider in detail LDA+DMFT(QMC) results for spectral function A(k,u) for 
SrVOs obtained in Ref. [92]. Owing to ideal cubic lattice symmetry self-energy matrix 
E(a;) is diagonal and all diagonal elements are the same for all t2 g orbitals. Spectral 
function is defined by imaginary part of Green function ImG(k, u), i.e. in fact by self- 
energy on real axis. This self-energy was calculated by numerical solution of Dyson's 
equation for the known interacting and bare Green functions as described in the Appendix 
of Ref. [T22] . 

In the Fig. [TH] this self-energy £(o>) is plotted as a function of real frequencies. It 
is essentially asymmetric with respect to the Fermi level, as could be assumed from the 
asymmetry of LDA density of states and band filling 1/6. At energies to ~ ±1.5 eV the 
real part of self-energy has extrema, corresponding to transition region from quasiparticle 
peak to lower and upper Hubbard bands. Two extrema in the imaginary part of the 
self-energy coinciding with ReS zeros determine energy positions of lower and upper 
Hubbard bands (see Fig. [T8|) . 

Asymmetric quasiparticle peak in the density of states (DOS) is situated in the energy 
range -0.8-1.4 eV Fig. [18j We see that imaginary part of self-energy ImS(w) is sufficiently 
small for these energies, while the real part can be roughly approximated with the dashed 
straight line shown in Fig. [T9]) . The slope of this line defines mass renormalization value 

7 Here we remind that real part of self-energy is connected with its imaginary part via Kramers-Kronig 
relation [5T] . 
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Figure 19: Real (black solid line) and imaginary (grey line) parts of self-energy £(o>) 
obtained from LDA+DMFT(QMC) calculation for V3d-t 2g orbitals of SrV0 3 [92]. In the 
inset self-energy in the vicinity of the Fermi level is shown in more details. Dashed lines 
give the slope of ReE(a;) far away and close to the Fermi level. 

Z = m* /m = 1 — dRc ^j^ \ uj=o — 1-9. Such Z value agrees with one obtained from QMC 
data in Matsubara frequencies: m* /m = 1 — Im ^°- > w 2, where Uq is "zeroth" Matsubara 
frequency. This value of renormalization is in accord with the value m* jm = 2.2 from the 
works [1371 138J as well as with experimental estimate from ARPES data |139j . 

In the inset in Fig. [19] it is visible that Fermi liquid behavior of self-energy ImS(oj) ~ 
— u 2 , together with Re£(o;) ~ —u, is fulfilled only in the interval from -0.2 to 0.15 eV. 
The slope of Re£(a;) in immediate proximity to the Fermi level is steeper than in wider 
energy interval (dashed line in Fig. IT9l) . Thus Fermi liquid mass renormalization value is 
larger than m* /m = 1.9 and is equal to iriy owE /m = 3 (dashed line on the inset of Fig. [T9l. 
On the edges of Fermi liquid regime, sharp bends of ReS(o;) at energies u = ±0.25 eV are 
seen. As the border of Fermi liquid regime we can consider energies, where the self-energy 
behavior starts to differ from Im£(u;) ~ — u 2 , which because of Kramers-Kronig relation 
corresponds to Re£(u;) ~ —u. Deviation from the square behavior of ImE at energies of 
the order of u — ±0.25 eV immediately leads to cusps in the Re£(a;). 

If the self-energy on the real axis is known one can compute spectral function A(k, u) 
and also the quasiparticle dispersion determined by momentum dependence of its maxima. 
In Fig. [20] we show the map of the spectral density for SrVC>3 obtained in Ref.[92]. In 
this multiband system (with degenerate bands) further analysis is similar to that for 
the single band case of section 13.3.21 White dots denote dispersion curves E nk obtained 
from LDA+DMFT calculation for SrVC>3. In the narrow vicinity of the Fermi level they 
coincide with LDA band structure e„k (line 1) renormalized by Fermi liquid factor Zfl = 
0.35, so that En^ = Z^e^ (line 2). Outside the Fermi liquid region dispersion curves 
correspond to LDA band structure with different renormalization factor: E n ^ = Zcp€ n k + 
c± (line 3), where Zqp = 0.64, c + =0.086 eV, c_ =0.13 eV. Along high-symmetry T-M 
and T-R directions in the Brillouin zone, transition between these two regimes leads to 
formation of kinks in the effective dispersions at energies w* i+ =0.22 eV and aj* 5 _=-0.24 
eV. These kinks are marked with arrows on the right side of Fig. [20] which corresponds 
to the area surrounded with the white rectangle on the main part of the figure. On the 
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Figure 20: Kinks in the dispersion relation E n y_ (white dots), for SrVC>3 obtained from 
LDA+DMFT calculation. Line 1 - LDA band structure e n t; line 2 - LDA bands renor- 
malized with Fermi liquid mass renormalization factor, lines 3 - band structure within the 
intermediate regime. Kinks are marked with arrows on the right panel, which corresponds 
to the area restricted with white rectangle on the main part of the figure. 

contour plot of spectral function A(k, uj) it is seen that in the energy region sufficiently far 
away from the Fermi level, spectral function keeps explicit k-dependence, despite pretty 
large damping value, replacing traditional band structure picture for systems with strong 
electron-electron correlations. 

Kinks of electronic nature were discovered in this system in ARPES experiments |139] 
in T-M direction at energies of the order of 0.25 eV, which agrees quite well with results 
of LDA+DMFT(QMC) calculations. 

4.4 Electronic structure of copper oxides in the pseudogap state: 
LDA+DMFT+S. 

Pseudogap state as was already pointed above is one of main anomalies of the normal 
state of high-T c cuprates and it is thought that clarification of its physical nature is the 
key point to understand high-temperature superconductivity mechanism [381 EHJ 1140] . 
Most powerful tool to investigate this state in recent years became angular resolved pho- 
toemission (ARPES). During last ten years in this area there was a remarkable progress 
related to significant growth of ARPES resolution both in energy and in momentum space 
|141[ 1142] . From ARPES data the Fermi surface (FS) shape, quasiparticle dispersion and 
damping, even self-energy can be directly restored |141[ 1142] . This allowed to study in 
detail formation of the pseudogap, "shadow" bands, quite unusual phenomena of Fermi 
arcs formation, interlayer hybridization effects (bilayer splitting) in double layer systems, 
|141[ 1142] , to determine qualitative distinctions between electron and hole doped cuprates 
[1411 1142] . The purpose of theory is an explanation of all these peculiarities and this 
problem is much complicated by rather strong electronic correlations, typical for these 
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systems and making doubtful the standard band theory and Fermi liquid approach. 

In this section we shall demonstrate that an account of AFM short range order fluctu- 
ations is in principle enough to describe a number of ARPES experiments in real systems. 
To this purpose we use LDA+DMFT+E hybrid computational scheme [143|. 1144} 1145} 
I146[ 11471 1148] . On the one hand this scheme inherits all advantages of LDA+DMFT 
[EH [__, US] HHJ EO] , i-e. the combination of single electron first principle density functional 
theory within the local density approximation (DFT/LDA) |108[ 1149] with dynamical 
mean- field theory (DMFT) for strongly correlated electrons [9j [3j HJ [7]. On another 
hand this scheme allows to consider nonlocal correlations by introduction of momentum 
dependent self-energy, while the usual self-consistent set of DMFT equations is preserved 
[231 124"! [25] . To solve effective single impurity problem of DMFT in the works described 
below we employed numerical renormalization group (NRG) [T4"| [15] . 

Such computational scheme fits very well to describe electronic properties of high-T c 
cuprates in normal (underdoped) state. Firstly, all material specific model parameters 
of physically relevant Cu-3d x 2 — y 2 orbital can be obtained from LDA calculations. 
Secondly, stoichiometric cuprates are antiferromagnetic Mott insulators with U ^> W 
(U - local Coulomb interaction W - conduction band width), so that correlation effects 
there are very important. At finite doping (at least up to optimal doping) cuprates are 
typical strongly correlated metals and DMFT stage of the computational scheme allows 
one to account for strong electronic correlations. Finally, to study "antiferromagnetic 
scenario" of pseudogap formation we introduce into the standard LDA+DMFT scheme 
k- dependent self-energy Ek, describing nonlocal correlations induced by (quasi) static 
Heisenberg spin fluctuations of short range AFM order [HJ H2] . 

In the framework of LDA+DMFT+E approach we performed calculations for a series 
of high temperature superconductors: hole doped Bi 2 Sr 2 CaCu 2 8 _ ( 5 (Bi2212) [143J and 
La 2 _ :r Sr :E Cu04 (LSCO) [HI], and also for electron doped Nd2_a.Cea.CuO4 (NCCO) [Hg 
H_B] and Pr 2 _ a; Ce :r Cu04 (PCCO) jl_7) . LDA+DMFT+E calculation results for Fermi 
surfaces and spectral functions can be compared with ARPES data for quasiparticles 
bands and experimental Fermi surface maps. 

Crystal structure of EU2212 [143], NCCO [1451 EE] and PCCO [147] have tetrago- 
nal symmetry with the space group I4/mmm, while LSCO has distorted orthorhombic 
structure Bmab [144j . In more details crystallographic data used in LDA+DMFT+E 
calculations are presented in Refs. P331 Ell E3 ESI E3 EE]- 

It is well known that physical properties of cuprates in many respects are determined 
by quasi-two-dimensionality of their electronic properties. From this point of view, the 
main interest is attracted to electronic states of Cu0 2 plane, where we are dealing with 
partially filled antibonding Cu-3d(x 2 —y 2 ) orbitals with dispersion crossing the Fermi level. 
In tight-binding approximation this dispersion has the following form: 

e(k) = —2t (cos k x a + cos k y a) — At' cos k x a cos k y a (66) 
— 2t" (cos 2k x a + cos 2k y a) — 2t'"(cos k x a cos 2k y a + cos 2k y a cos k y a). 

Here t, t', t", t'" are Cu-Cu transfer integrals in first four coordination spheres in the Cu0 2 
plane, a is the lattice constant. Values of these effective transfer integrals calculated with 
the use of Wannier functions obtained within the N-th order muffin-tin orbitals method 
(NMTO) of Ref. [m] are listed in the Table HJ In the following we shall exploit LDA 
calculated effective antibonding Cu-3d(x 2 —y 2 ) band as a "bare" one in LDA+DMFT+E 
calculations. 
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Table 1: Calculated model energy parameters (eV) and experimental correlation length 
£. First four Cu-Cu transfer integrals in the Cu02 plane t, t', t", t'"; effective interlayer 
transfer integral t±, local Coulomb interaction U and pseudogap amplitude A. 





t 


t' 


t" 


t'" 


t± 


U 


A 




Bi2212 


-0.627 


0.133 


0.061 


-0.015 


0.083 


1.51 


0.21 


10a 


NCCO 


-0.44 


0.153 


0.063 


-0.01 




1.1 


0.36 


50a 


PCCO 


-0.438 


0.156 


0.098 






1.1 


0.275 


50a 


LSCO 


-0.476 


0.077 


-0.025 


-0.015 




1.1 


0.21 


10a 



In double layer systems, e.g. in Bi2212, hopping between two neighboring planes is 
also important. In tight-binding approximation an expression for corresponding interlayer 
dispersion derived in Ref. |150] has the form: 



t_i_(k) = -^(cosA^a — cos k y a) 2 (67) 

The value of t± is given in Table [TJ Consideration of interlayer hopping and "bilayer 
splitting" effects requires certain generalization of LDA+DMFT+E computational scheme 

To perform DMFT calculations one should also calculate a value of onsite Coulomb 
interaction. The value of this interaction for effective Cu-3d(x 2 — y 2 ) orbital obtained 
within the constrained LDA method |117] is also given in Table [U 

To account for AFM spin fluctuations we employed two-dimensional model of the pseu- 
dogap state [Ml 02], generalized for DMFT+E calculations [23J [25] . Additional "external" 
k-dependent self-energy Ek [231 [25] describes nonlocal correlations caused by (quasi)staticQ 
AFM spin fluctuations. 

To specify Ek it is necessary to know two important parameters - the pseudogap 
amplitude A, giving energy scale of fluctuating SDW and correlation length £. The 
value of A were calculated as described in Refs. [23J, [251 1 143] . The values of correlation 
length were taken in accordance with values obtained in neutron scattering experiments 
for NCCO [151] and LSCO [152] . The values of A and £ used for all systems under 
consideration are also listed in Table [U To solve effective Anderson single impurity 
problem in DMFT we used numerical renormalization group (NRG [THUS]). Temperature 
in DMFT(NRG) calculations was chosen to be 0.011 eV and electron or hole concentration 
(doping level) was taken to be 15%. 

LDA+DMFT+E calculations produce a clear picture of "hot spots" behavior in the 
spectral function and on maps of the Fermi surfaces for electron doped systems [1451 11461 
I147j , while for hole doped systems only Fermi arcs arise [143[ 1144] . 

In Fig. |2~T1 LDA+DMFT+E we show spectral functions along 1/8 part of bare Fermi 
surface from nodal point on the diagonal of the Brillouin zone (upper curve) downto 
antinodal point at the boundary of the zone (lower curve). Results for Bi2212 are shown 
on the left panel and for NCCO on the right panel of Fig. [21] For both systems in 
nodal direction quasiparticles are well defined — sharp peak of spectral function situated 



8 Quasistatic approximation for AFM fluctuations necessarily restricts this approach to rather high 
temperatures (and energies not very close to the Fermi level) [41] [42] . Thus we can not judge about the 
nature of low temperature (low energy) damping which is determined by dynamical (inelastic) scattering 
processes. 
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Figure 21: LDA+DMFT+S spectral functions for Bi2212 (leftpanel) and NCCO (right panel) 
along the "bare" Fermi surface in the l/8th of the Brillouin zone. Black line corresponds to the 
"hot spot" [US [T36] 



practically on the Fermi level is clearly seen. As one moves to antinodal point quasiparticle 
damping grows reaching the maximum at the "hot-spot" and the peak of spectral density 
moves away from the Fermi level. This behavior is in complete agreement with results of 
Refs. |153l 1154] (comparison with experiment see in [1451 1146"] ). From LDA+DMFT+E 
results shown in Fig. [21] it is directly seen, that for Bi2212 antinodal states are formed by 
low energy edge of the pseudogarj^, while for NCCO by high energy edge. For Bi2212 we 
also observe bilayer splitting of quasiparticle peak which is related to the presence of two 
Cu02 planes in the elementary cell. 

The "hot spots" for NCCO are located closer to the Brillouin zone diagonal |145[ 
1146] . This can be seen from black lines on Fig. [21] which correspond to the "hot spots" . 
Moreover, correlation length in NCCO is much larger than in Bi2212. Thus for NCCO (in 
contrast to Bi2212) in the antinodal direction quasiparticles again are rather well defined. 
For Bi2212 scattering near by Brillouin zone boundaries is strong everywhere and instead 
of "hot spots" picture we observe quite strong "destruction" of the Fermi surface close to 
these boundaries. Qualitatively the same picture is observed also in LSCO. 

In Fig. [22] on the upper panel, LDA+DMFT+ZI Fermi surface maps in the quarter 
of the Brilloun zone for Bi2212 (on the left) and NCCO (on the right) are presented. 
In Bi2212 we observe strong "destruction" of Fermi surface by scattering by pseudogap 
fluctuations close to Brillouin zone boundaries^!. On the contrary, in NCCO the Fermi 
surface is almost recovered close to the Brillouin zone boundaries. On the other hand 
Fermi arc around nodal direction in Bi2212 is clearly pronounced, while for NCCO it is 
noticeably smeared. This is another consequence of the fact that in NCCO "hot spots" 

9 Especially clear it is visible for the case of smaller correlation length £ = 5a considered in Ref. [143] . 
10 Analogous behavior in entire accordance with ARPES results is realized also in another hole doped 
system - LSCO. The LDA+DMFT+S calculations for this system were performed in [144) 
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Figure 22: LDA+DMFT+S Fermi surface maps obtained in Refs. [HSlfHBj for Bi2212 (upper 
left panel) and NCCO (upper right panel) in the quarter of the Brillouin zone, (&£, ky are given 
in units of ir/a). Experimental Fermi surface for Bi2212 (lower left panel |155| ) and NCCO 
(lower right panel [153] ) . 



are located closer to the Brillouin zone diagonal. A bit larger value of the pseudogap 
amplitude A also favours the stronger damping of Fermi arcs in NCCO. One should note 
the appearance of the "shadow" Fermi surface, which is much more intensive for NCCO. 

Qualitatively the same Fermi surfaces were observed experimentally in real Bi |155] 
and Nd [153J systems (lower panel of Fig. [221) . Thus, the distinction of Fermi surface 
maps for these systems is related mainly to the distinction of band structure parameters 
of these materials. In particular, LDA Fermi surfaces of NCCO are more bended and 
"hot spots" appear to be rather far from the Brillouin zone boundaries, consequently the 
Fermi surface in the vicinity of these boundaries almost does not feel scattering by AFM 
fluctuations. In Bi2212 LDA Fermi surface is rather close to the Brillouin zone boundaries 
and (7r/a,0) point, so that "hot spots" are also close to this point. Thus in Bi2212 they 
are more "washed away" by strong pseudogap scattering close to (n/a,0) point and are 
not observed. "Hot spots" in NCCO are more vivid also because of much larger value of 
correlation length of fluctuations. 

Not less graphic results were obtained in LDA+DMFT+E calculations and ARPES 
experiments for [147] . In Fig. [23] we show PCCO Fermi surface map (panel (a) - 
LDA+DMFT+S results, panel (b) - experimental ARPES data). Fermi surface here 
is clearly distinguishable only near Brillouin zone boundaries and around (n/a/2,7r/a/2) 
point (Fermi arc). Again, as in NCCO we observe "destruction" of the Fermi surface in 
"hot spots" , located at the intersection of Fermi surface and its AFM shadow "replica" is 
detected. This " destruction" of the Fermi surface is due to the strong electron scattering 
by AFM spin (pseudogap) fluctuations. The "shadow" Fermi surface is observed, as it 
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Figure 23: Fermi surface map for PCCO. (a) — results of LDA+DMFT+E calculation. White 
rectangle on the panel (a) shows the part of reciprocal space where ARPES measurements (panel 
b) were done. Lower left corner coincides with X-point (tt/cl, 0) of the Brillouin zone |147j . 

happens in the case of AFM doubling of the lattice period. However, since there is no 
long-range order in the underdoped region, in which we are interested, this "shadow" 
Fermi surface is strongly eroded. Fermi surface of PCCO is very similar to that ob- 
served in Nd2_ 2: Ce 2: Cu04 (NCCO), which belongs to the same family of superconduc- 
tors p^miT53] . 

Let us compare (see Fig. I24"j) theoretical (upper panel) and experimental (lower panel) 
quasiparticle dispersions along most characteristic cuts of the Fermi surface shown in 
Fig. [23J Theoretical data are multiplied by the Fermi function with a temperature of 30K 
and convoluted (in energy) with a Gaussian distribution to simulate the experimental 
resolution. Cut 1 crosses the quasiparticle and the "shadow" Fermi surface near the Bril- 
louin zone boundary. Correspondingly, here it is possible to detect "fork" -like structure, 
formed by suppressed "shadow" band and much better defined quasiparticle band. This 
structure corresponds to the beginning of formation of the Fermi surface cylinder around 
( pi/ a, 0) point. Cut 2 passes exactly through the " hot spot". Here we see a strong 
suppression of the quasiparticle band near the Fermi level. Cut 3 crosses the Fermi arc 
and we can see fairly well-defined quasiparticle band. However, the " shadow" band of low 
intensity is also present. In the case of long-range AFM order and a full doubling of the 
period, the Fermi surface and its "shadow" form a closed "pocket" , of the Fermi surface 
around the (n/2a, 7r/2a) point, while in the present case a part of the pocket formed by 
"shadow" band is strongly blurred. One can see that there is a good agreement between 
calculated and experimental data. 

As was already noted, within LDA+DMFT+E scheme two-particle properties can 
also be calculated |26j, which allowed to investigate optical conductivity of Bi and Nd 
cuprates |145[ 1146] , also demonstrating significant differences in the effects of pseudogap 
fluctuations. In particular, in optical conductivity of NCCO, in qualitative agreement with 
experiment |157j , we observe a characteristic pseudogap dip and a smooth maximum due 
absorption through the pseudogap at frequencies ~ 2 A. However, in optical conductivity 
of Bi2212 characteristic pseudogap structure practically does not occur neither in theory 
or in experiments |158j . which is related to sufficiently small values of A and fluctuation 
correlation length in this system. 

Let us summarize some of our conclusions. For all the systems studied, LDA+DMFT+S 
calculations show, that Fermi-liquid behavior persists only rather far away from the "hot 
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Figure 24: Momentum - energy distribution curves for characteristic cuts of the Brillouin zone 
plotted in Fig. [23] (upper panels — theoretical data, lower panels — experimental photoemission 
intensity). For the cut 1 we show momentum distribution curve (MDC) integrated over energy 
window of the width of 60 meV around the Fermi level. Analogous MDC for the cut 2 (through 
the "hot spot") shows ARPES intensity suppression as compared with MDC for the cut 3, which 
is located quite far away from the "hot spot" [145} 1146] , 

spots" (nodal direction), and "destruction" of the Fermi surface is observed near the "hot 
spots" . This destruction is due to strong scattering of correlated electrons on short-range 
order AFM (pseudogap) fluctuations. Comparison of ARPES data and LDA+DMFT+S 
calculations shows the existence of quite distinct "hot spots" in the behavior of the spec- 
tral density and maps of the Fermi surface in electron doped systems, in contrast to hole 
doped systems, where we only observe only a strong "destruction" of the Fermi surface 
near Brillouin zone boundaries and Fermi arc around its diagonal. There are several rea- 
sons for this difference: (1) "hot spots" in electron doped systems are located closer to 
the center of the Brillouin zone, (2) correlation length of AFM fluctuations in electron 
doped systems is bigger; (3) the width of the pseudogap in the electron doped systems is 
also larger than in hole doped ones. Experimental and theoretical results discussed here 
clearly confirm the AFM scenario of the formation of the pseudogap in both hole-doped 
[T4"31 QH], and electron doped cuprates [145] HM\ H47] . 

5 Conclusion 

In this review we discussed DMFT + E generalization of the standard dynamical mean- 
field theory (DMFT), which allows to include non-local correlations or additional (rel- 
ative to Hubbard one) interactions (in principle of any type), while remaining within 
the single-impurity picture of DMFT and retaining the same set of self-consistent DMFT 
equations. The basic approximation of this method is the neglect of interference contribu- 
tions of DMFT diagrams and additional interactions included into the analysis. Precisely 
this (strictly speaking not completely controllable) approximation allows to preserve the 
overall structure of DMFT equations, which permits to solve DMFT+E equations with 
well developed methods used in the standard DMFT. It must be emphasized that the 
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self-consistent account of additional interactions at every step of DMFT loop leads to a 
rather complicated procedure, equivalent to the summation of infinite classes of diagrams. 

The proposed approach proved to be versatile enough to be applied to a number 
of problems in systems with strong electron correlations - from semi-phenomenological 
account of non-local short range order pseudogap fluctuations to the self-consistent scheme 
for metal-insulator transition in the disordered Hubbard-Anderson model and account 
of the effects of electron-phonon interaction in electronic spectra of strongly correlated 
systems. A remarkable feature of DMFT+S approach is the possibility to study, along 
with one-particle characteristics, also the two-particle properties, i.e., in principle, any 
kind of response functions (optical conductivity, magnetic susceptibility, charge screening, 
etc.). The universality of the method allows one to hope for its successful application in 
a number of future problems. 

Discussing all problems under consideration, one should keep in mind that in many 
respects similar physical results can be obtained with more sophisticated approaches, using 
these or other methods of direct numerical simulation. For example, similar results for the 
formation of the pseudogap in the single-particle characteristics of the two-dimensional 
Hubbard model were obtained in the cluster generalizations of DMFT [2T| [22] . However, 
these methods have specific limitations (e.g. of cluster size) and are still not widely used to 
calculate the two-particle properties, such as the general response functions, in particular, 
the optical conductivity. DMFT + X approach has obvious advantages, associated with 
savings of computational resources. It requires a significantly lower cost of computational 
time, and its advantage in calculating the two-particle response functions is quite obvious. 
This opens up additional opportunities for the systematic comparison of various types 
of non-local fluctuations or additional interactions and their influence on the electronic 
properties of strongly correlated systems, providing intuitively clear path to the analysis 
of experiments and theoretical results obtained via more complicated schemes. 

Rather simple generalization of our computational scheme enabled us to formulate 
also the generalized LDA+DMFT+S approach, that allows to perform calculations of 
all of the effects discussed above for real compounds of transition elements with strong 
electronic correlations. One can expect that these calculations will be useful in analyzing 
and explaining the new experimental data. 
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